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Preface 


My  quest  for  a  thesis  topic  began  with  an  interest  in  detecting  space  objects  us¬ 
ing  a  space-based  sensor.  My  advisor,  Dr.  Michael  C.  Roggemann,  worked  with  me 
on  defining  the  scope  of  my  research.  During  a  visit  to  Phillips  Lab,  Dr.  Roggemann 
discovered  the  Lasers  and  Ima^ng  branch  (PL/LIMI)  was  interested  in  hyperspectral 
imaging  of  space  objects  for  Space  Object  Identification  applications.  Hyperspectral 
imaging  involves  viewing  an  object  in  very  narrow  wavelength  bandwidths.  The  in¬ 
formation  stored  in  each  wavelength  bandwidth  forms  a  spectral  signature,  which  can 
possibly  be  used  to  identify  an  object.  My  thesis  evolved  into  hyperspectr<d  analysis 
of  space  objects,  and  addressed  the  question  if  enough  target  signal  was  present  in 
the  narrow  wavelength  bandwidths  to  conduct  hyperspectral  measurements.  I  owe  a 
great  deal  of  thanks  to  Dr.  Roggemann  for  his  guidance,  insight,  patience,  and  sense 
of  humor,  I  can  honestly  say  I  enjoyed  my  research.  I  also  thank  Dr.  Theodore  E. 
Luke,  my  co-advisor,  for  asking  difficult  questions  that  probed  my  understanding. 
To  my  wife  Danielle  and  son  Peter,  I  love  you  very  much,  and  thanks  for  everything. 
Finally,  I  thank  God,  who  makes  cdl  things  possible. 

I  dedicate  this  thesis  to  General  James  H.  Doolittle,  pioneer  aviator,  1923  grad¬ 
uate  of  the  Air  Service  Engineering  School  (AFIT’s  predecessor),  and  leader  of  the 
Tokyo  raid.  General  Doolittle  advocated  effective  use  of  air  power  for  the  defense  of 
the  United  States.  I  pray  I  have  the  same  vision  for  space  power  as  General  Doolittle 
had  for  air  power. 


General  James  H.  Doolittle,  1896-1993 


Daniel  Hrovat 
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Abstract 

Hyperspectral  SLnalysis  of  space  objects  is  being  considered  as  a  me2uis  of  en¬ 
hancing  Space  Object  Identification  (SOI)  capabilities.  This  study  develops  a  math¬ 
ematical  simulation  based  on  fundamental  radiometry  principles  to  evaluate  the  fea¬ 
sibility  of  hyperspectral  analysis  of  space  objects.  The  study  uses  Signal  to  Noise 
Ratio  (SNR)  as  the  basis  for  the  feasibility  evaluation,  and  addresses  the  question: 
Is  there  enough  target  signal  in  narrow  wavelength  bandwidths  to  conduct  hyper- 
spectr2d  measurements  of  space  objects?  The  analysis  characterizes  the  target  signal 
and  noise  sources  using  fundamental  raidiometry,  and  develops  a  flexible  radiometry 
model  which  calculates  the  target  signal  count  rate,  and  SNR  for  a  specified  target¬ 
sensor  configuration.  The  model  was  applied  to  four  case  studies,  and  revealed:  1) 
Enough  tary**'  ignal  is  present  in  narrow  wavelength  bandwidths  for  hyperspectral 
measurements;  2)  Differences  exist  in  signal  strength  as  a  funcion  of  wavelength  and 

Object  Identification  applications, 
further  investigation. 


bandwidth  which  may  serve  as  a  basis  for  Space 
Hyperspectral  emalysis  of  space  objects  warrants 
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HYPERSPECTRAL  ANALYSIS  OF  SPACE  OBJECTS; 
SIGNAL  TO  NOISE  EVALUATION 


L  Introduction 

Hyperspectral  analysis  of  space  '  Ejects  is  being  considered  as  a  means  of  en¬ 
hancing  Space  Object  Identification  (SOI)  capability.  This  study  examines  the  feasi¬ 
bility  of  hyperspectral  analysis  from  a  Signal  to  Noise  perspective.  This  introductory 
chapter  provides; 

•  Background 

•  Problem  Statement 

•  Approach 

•  Scope 

•  Summary  of  key  results 

•  Organization  of  the  thesis 

1.1  Background 

Hyperspectral  analysis  involves  viewing  an  object  in  narrow  bandwidths  of  the 
electromagnetic  spectrum.  The  wavelength  bandwidth  is  typically  on  the  order  of 
ten  nanometers  (4:79).  In  contrast,  multispectral  analysis  typically  uses  wavelength 
bandwidths  of  forty  nanometers  (5:456).  The  use  of  hyperspectral  analysis  to  iden¬ 
tify  a  space  object  requires  a  detailed  knowledge  of  the  spectral  signatures  of  the 
materials  used  to  make  the  object  (4:76).  A  report,  Matter.dat  Surface  Properties 
Data  Compilation.,  contains  the  spectral  reflectances  of  materials  commonly  us>  d  in 
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spacecraft  construction  (14).  Introduction  to  the  Physics  and  Techniques  of  Remote 
Sensing  describes  the  ideal  case  for  identifying  a  material  using  spectreil  analysis: 

...  if  a  certain  material,  or  family  of  materials,  is  the  only  one  which 
has  a  certain  spectral  feature,  such  as  2in  absorption  line  at  a  certain 
wavelength,  the  identification  becomes  simple.  The  identification  feature 
could  be  a  single  absorption  line,  or  an  association  of  lines.  (4:77) 

In  a  similar  way,  hyperspectral  analysis  of  space  objects  may  provide  a  specific  signa¬ 
ture  through  “an  association  of  lines”  unique  to  that  space  object.  This  hyperspec¬ 
tral  signature  would  enhance  Space  Object  Identification.  In  addition,  hyperspectral 
cuicilysis  may  provide  insight  into  a  space  object’s  material  composition.  The  two 
general  categories  of  hyperspectral  analysis  are  hyperspectral  imaging  zmd  hyper¬ 
spectrometry.  Differentiating  imaging  from  spectrometry  is  discussed  in  Chapter  II. 
Hyperspectral  analysis  is  being  considered  for  ground-based  eind  space-based  sensors 
(12)  (11).  The  ground-based  sensor  is  a  1.6  m  diameter  telescope  located  in  Maui, 
Hawaii.  A  proposed  space-based  sensor,  part  of  the  Space  Surveillance  Initiative 
(SSI)  research  program  (11),  may  have  a  hyperspectral  sensor  as  p2irt  of  its  mission 
suite. 

1.2  Problem  Statement 

Develop  a  mathematical  simulation  based  on  fundamental  radiometry  princi¬ 
ples  to  eveduate  the  feasibility  of  hyperspectral  analysis  of  space  objects. 

1.3  Approach 

The  key  metric  in  determining  a  sensor’s  feasibility  and  performance  is  the 
Signal  to  Noise  Ratio  (SNR).  The  SNR  is  a  function  of  the  target  signal,  noise, 
and  integration  time.  This  SNR  evaluation  addresses  the  question:  Is  there  enough 
target  signal  in  narrow  wavelength  bandwidths  to  conduct  hyperspectral  me2isure- 
ments?  The  approach  of  this  thesis  is  to  characterize  the  target  signal  and  noise 
sources  using  fundamental  radiometry.  The  integration  time  is  case  dependent,  and 
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is  approached  from  a  best  case,  worst  case  standpoint  to  bound  the  integration  time. 
The  equations  derived  from  the  radiometry  are  evaluated  in  a  mathematical  simula¬ 
tion  using  Mathcad  3.1  (6).  This  radiometry  model  is  applied  to  four  case  studies: 
two  involve  hyperspectral  imaging,  and  two  involve  hyperspectrometry. 

1.4  Scope 

1.  This  study  considers  the  radiometry  from  point  to  point,  no  atmosphere.  The 
radiometry  model  developed  for  the  SNR  evaluation  can  be  applied  to  ground- 
based  (first  order  approximation  with  no  atmosphere)  or  space-based  sensors. 

2.  The  sensor  is  tracking  the  target. 

3.  Background  noise  is  small  compared  to  the  target  signal,  so  background  noise 
is  negligible. 

4.  Spectral  reflectance  of  materials  is  not  linked  to  spatial  coordinates. 

5.  Direct  solar  irradiance  is  the  sole  source  of  energy  used  in  the  radiometry 
development. 

6.  The  detector  for  this  study  is  a  Chjurge  Coupled  Device  (CCD). 

1.5  Summary  of  Key  Results 

For  Case  (1)  and  Case  (2),  the  target  was  modeled  as  rocket  body  6  m  long 
and  2  m  in  diameter.  The  range  to  the  target  was  assumed  to  be  800  km.  For 
Case(l),  the  rocket  body  had  an  aluminum  surface.  For  Case  (2),  the  rocket  body 
was  painted  white.  For  Case  (3)  and  Case  (4),  the  targets  were  modeled  as  cylin¬ 
drical  communications  satellites  in  geosynchronous  orbit.  Range  to  the  targets  Weis 
36,000  km.  For  Case  (3),  the  satellite  was  7.6  m  long  and  3.64  m  in  diameter,  and 
has  an  antenna  2.4  m  in  diameter  on  top  of  the  cylinder.  For  Case  (4),  the  satellite 
is  4.77  m  long  and  2.2  m  in  diameter,  and  has  an  antenna  1.83  m  in  diameter  on  top 
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of  the  cylinder.  The  sensor  viewing  the  targets  for  all  cases  was  a  1.6  m  diameter 
telescope.  For  the  cases  studied: 

1.  Enough  signal  is  present  for  hyperspectral  imaging  and  hyperspectrometry  for 
narrow  wavelength  bandwidths  of  AA  =  20 nm,  10 nm,  and  5nm. 

2.  Analysis  of  the  radiometry  model  for  different  materials  and  objects  shows 
differences  in  signal  strength  as  a  function  of  wavelength  which  may  serve  as  a 
basis  for  Space  Object  Identification. 

1.6  Organization 

•  Chapter  I,  Introduction. 

•  Chapter  II,  Theory.  Establishes  the  theoretical  basis  for  the  radiometry  model. 

•  Chapter  III,  Results.  Outlines  the  procedures  for  running  the  radiometry 
model,  and  presents  the  results  from  the  four  case  studies. 

•  Chapter  IV,  Conclusions  t  d  Recommendations.  Provides  a  summary  of  the 
study,  the  key  results,  corollary  issues,  and  recommendations  for  further  re¬ 
search. 

•  Appendix  A.  Contains  the  Mathcad  code  used  to  develop  the  radiometry  model. 

•  Appendices  B-E.  Contains  the  graphical  output  for  the  four  case  studies. 

•  Appendix  F.  Contains  plots  of  spectral  reflectance  versus  wavelength  for  the 
materials  used  in  the  case  studies. 
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II.  Theory 

2. 1  Introduction 

The  main  metric  for  determining  the  feasibility  of  an  imaging  system  is  Signal 
to  Noise  Ratio  (SNR).  SNR  is  a  measure  of  a  system’s  ability  to  discriminate  between 
the  target  signal  and  noise.  Noise  originates  from  photon  shot  noise,  the  background, 
and  the  detector.  The  following  sections  address  the  theory  and  methodology  for 
determining: 

•  target  signal 

•  background  noise 

•  detector  noise 

•  integration  time 

•  SNR 

2.2  Target  Signal 

2.2.1  Introduction.  Radiometry  determines  the  target  signal  incident  on 
the  detector.  For  this  evaluation,  the  radiometry  is  derived  from  first  principles.  The 
result  of  the  radiometry  derivation  is  an  expression  for  the  average  count  rate  from 
the  detector.  The  count  rate  is  proportional  to  the  strength  of  the  target  signal,  and 
can  be  expressed  in  terms  of  target,  sensor,  and  source  parameters.  The  final  coimt 
rate  equation  includes  simplifications  used  in  the  Mathcad  radiometry  model. 

2.2.2  Radiometry.  Radiometry  determines  the  number  of  photons  incident 
on  the  detector.  The  detector  generates  a  count  rate  based  on  the  strength  of  the 
target  signal.  The  count  rate  from  the  detector(lO)  is 

‘t  = 
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where, 


=  count  rate  from  detector  (1/s) 

=  power  incident  on  the  detector  {W) 
e(A)  =  photon  energy  (J) 

»/(A)  =  detector  quantum  eflBciency 

A  =  wavelength  (m). 

The  photon  energy  is 

£(A)  =  ^  (2.2) 


where, 

h  =  Planck’s  constant  =  6.626  x  10“^^  J  s 

c  =  speed  of  light,  vacuum  =  3  x  10®  m/s. 

The  power  incident  on  the  detector  in  terms  of  detecting  optic  area  and  target 
irradiance  (10)  is 


$e(A)  =  SjiET{X)To{X)  (2.3) 


where, 

Sji  =  area  of  detecting  optic  (m^) 

^j.(A)=  irradiance  at  the  sensor  due  to  the  target  {W/rri^) 
ro(A)  =  optical  transmission. 

The  target  irradiance  at  the  sensor  (10)  is 

Ej.(X)  =  (2.4) 


where. 
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Lj'{X)=  target  radiance  {Wfm^sr) 

Sj'  =  target  area  (m*) 

r  =  range  to  target  (m) 

$l  =  target-sensor  angle. 

The  target-sensor  angle  (and  source-target  angle  described  below)  are  measured 
with  respect  to  the  surface  normal  of  the  target  (Figure  2.1). 


Figure  2.1  Target-sensor  and  source-target  angles 

To  express  the  radiance  of  the  target  in  terms  of  irradiance  from  a  source,  the 
target  is  assumed  to  be  Lambertian  (2:22),  and 

.  Iy(A)  =  MliR  (2.5) 

TT 

where, 

Mj'{\y=  target  exitance  (IF/m^). 

The  target  exitance  (10)  equals 


Mj'{X)  —  EsourceiX)  p{X) 


(2.6) 


where, 


EsaurceW  =  irradiance  of  the  source  at  the  target  {Wfm'^) 
p{X)  =  spectral  reflectivity. 

For  this  analysis,  the  sun  is  the  only  source  of  irradiamce  on  the  target 

Esource{X)  —  Esun{X)  cos(02)  (2-7) 


where. 


Esun{X) 

h 


=  irradiance  of  the  sun  at  the  target  (VF/m^) 
=  source-target  angle 


The  resulting  expression  for  the  count  rate  is 

=  ?Zcas{ei)cos(02)ST  l^i^pWroWvWXdX  (2.8) 

Although  h,  c,  and  tt  are  not  dependent  on  wavelength,  retaining  these  constants 
in  the  integral  will  allow  Equation  2.8  to  be  calculated  in  terms  of  a  photon  flux, 
described  in  Section  2.2.3.  The  elements  of  Equation  2.8  can  be  categorized  accord¬ 
ing  to  target,  sensor  and  source  parameters.  (Table  2.1)  The  following  subsections 
describe  each  of  the  parameters  listed  in  Table  2.1. 


Taurget 

Sensor 

Source 

H 

AA 

ro(A) 

vW 

Esun{X) 

h 

Table  2.1  Target,  Sensor,  and  Source  Parameters 
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2.2.2. 1  Target  Parameters.  Target  Area  {Sj'):  As  target  area  in¬ 
creases,  the  count  rate  increases.  Target  size  can  very  from  a  few  centimeters  to 
greater  than  10  m  in  length. 

Spectral  Reflectivity  (p(A)):  Spectral  reflectivity  varies  for  each  material  with  0  < 
p(A)  <  1,  and  is  dependent  on  wavelength. 

Range  to  target  (r):  The  count  rate  decreases  by  the  square  of  the  range  to  the 
target. 

Target-sensor  angle  {$i  ):  For  this  analysis,  the  taxget-sensor  angle  has  a  range  of 
0  deg  <  B\  <90  deg,  therefore  1  >  cos{^l)  >  0  .  The  count  rate  is  a  maximum  for 
normal  incidence,  when  =  0  deg.  This  amadysis  only  considers  the  sensor  viewing 
one  side  of  the  target. 

2. 2. 2. 2  Sensor  Parameters.  Wavelength  (A):  The  wavelength  identi¬ 
fies  which  portion  of  the  electromagnetic  spectrum  the  sensor  uses.  Wavelength  for 
the  diflFerent  case  studies  varies  from  300  nanometers  (nm)  as  the  lowest  wavelength, 
to  850  nm  as  the  highest  (9).  The  choice  of  detector  will  limit  the  wavelengths  of  the 
sensor.  A  thinned  CCD  with  an  antireflective  coating  will  have  a  quantum  efficiency 
greater  than  50  percent  for  the  wavelength  region  of  300  nm  <  A  <  900  nm  (7). 
Wavelength  bandwidth  (A  A):  The  wavelength  b2uidwidth  gives  the  width  of  the 
divisions  within  the  wavelength  interval.  For  this  analysis,  wavelength  beindwidths 
support  Phillips  Lab’s  goal  of  hyperspectral  imaging  at  wavelength  bandwidths  of 
twenty,  ten  and  five  nanometers  (13). 

Area  of  detecting  optic  {Sj<):  As  the  area  of  detecting  optic  increases,  the  count  rate 
increases. 

Optical  transmission  (  to(A)):  Although  dependent  on  wavelength,  the  optical  trans¬ 
mission  will  not  vary  significantly.  The  optical  transmission  was  set  as  a  constcint. 
To  =  .98. 

^  Since  the  wavelength  detected  by  the  sensor  is  a  subset  of  the  wavelength  emitted  by  the  source, 
wavelength  is  categorized  as  a  sensor  parameter. 
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Detector  quantum  efficiency  (»?(A)):  The  quantum  efficiency  ranges  from  0  <  r}{X)  < 
1.  For  a  thinned  CCD  detector  with  an  antireflective  coating,  the  peak  quantum  ef¬ 
ficiency  is  approximately  70  percent  at  750  nm  (7).  For  the  wavelengths  used  in  this 
analysis,  the  quantum  efficiency  will  not  drop  below  50  percent.  Since  quantum  effi¬ 
ciency  between  the  maximum  and  minimum  only  differs  by  20  percent,  the  quantum 
efficiency  was  set  as  a  constant,  rj  =  .6. 


2.2.2.S  Source  Parameters.  Irradiance  of  the  sun  (F^sun(A)  ):  To 
determine  the  irradiance  of  the  sun  at  the  target,  the  sun  is  modeled  as  a  5770  K 
blackbody  (15).  The  spectral  exitance  from  the  sun  is  given  by  Planck’s  formula 
(4:35) 


where, 


MbbW 


2irh(^ 

^'hKxfer)-!, 


(2.9) 


h  =  Planck’s  const2mt  =  6.626  x  10“^'^  J  s 

c  =  speed  of  light,  vacuum  =  3  X  10®  m/s 

X  =  wavelength  (m) 

kg  =  Boltzmann’s  constant  =  1.38  x  10“^®  J/K 
T  =  temperature  {K) 

The  irradiance  cf  the  sun  at  the  target  is  given  by  the  inverse  square  law  (2:25) 


EsunW  = 

\  Tsun 


rsun  ) 


(2.10) 


where. 


MbbW  =  exitance  from  the  sun  (W/m^) 

Rsun  =  radius  of  the  sun  =  6.96  X  10®  m 
rsun  =  distance  to  the  sun  =  1.5  x  10^ 
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Source-target  angle  (^2)'  analysis,  the  source-target  angle  has  a  range  of 

0  deg  <  612  <  90  deg,  therefore  1  >  003(62)  >  0.  The  count  rate  is  a  maximum  for 
normal  incidence,  when  62  =  0  deg.  This  analysis  only  considers  the  irradiance  on 
one  side  of  the  target. 

2.2.3  Target  Count  Rate  Equation.  To  simplify  the  count  rate  equation 
(Equation  2.8),  the  optical  transmission  and  quantum  efficiency  will  not  depend  on 
wavelength,  so  the  form  of  the  count  rate  equation  is 

^  cos(6i)  003(62)  To  Sjr/  I  )  p{X)Xd\  (2.11) 

A  major  difficulty  with  this  equation  is  the  number  of  parameters  that  vary.  Tair- 
get  area,  range,  encounter  angles,  and  area  of  detecting  optic  vau-y  depending  on 
the  scenario.  A  means  of  simplifying  this  equation  is  through  a  hypernormalized 
parameter 

ft  =  003(62)  003(62)  To  Sj’ T?  (2-12) 

The  hypernormadized  parameter  models  a  wide  range  of  cases,  and  allows  the  model 
to  focus  on  the  outcome  —  the  tauget  count  rate.  Chapter  III  contaiins  a  section 
which  illustrates  the  usefulness  of  the  hypernormalized  parameter.  With  the  hyper- 
normalized  parameter,  the  count  rate  equation  simplifies  to 

!I  =  n  /^^f221p(A)AdA  (2.13) 

q  J  heir 

To  extract  the  spectrail  reflectivity  from  the  integral,  the  spectral  reflectivity  values 
from  a  given  material  are  treated  as  a  vector.  Each  element  of  the  vector  corresponds 
to  the  mean  wavelength  of  the  integral.  For  instance,  for  limits  of  integration  of 
445  nm  to  455  nm  (AA=  10  nm),  the  spectral  reflectivity  corresponds  to  wavelength 
450  nm.  With  the  spectral  reflectivity  removed  from  the  integral,  the  count  rate 
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equation  is 


9 


n  p 


J 


hen 


(2.14) 


The  remaining  terms  in  the  integral  are  the  irradieince  of  the  sun,  the  wavelength,  and 
constants  h,  c,  and  n.  From  Equation  2.10,  the  irradiance  of  the  sun  is  determined 
by  Planck’s  formula  and  the  inverse  square  law.  By  including  the  other  terms  in 
the  integral  (/i,c,A),  the  resulting  output  is  a  photon  flux,  P,  in  units  l/(m’ s).  The 
resulting  expression  for  the  target  count  rate  is 

—  =  Up-  (2.15) 

q  n 

where, 

n  =  hypernormalized  parameter  (radm^) 

p  =  spectral  reflectivity 

P  =  photon  flux  (l/(m^  s) 

2.3  Background  Noise 

Background  noise  sources  include  the  celestial  background,  earthlimb,  aurora, 
airglow,  and  zodiacal  light  (3:46).  These  noise  sources  will  decrease  SNR.  The  back¬ 
ground  signal  is  weak  compared  to  the  target  signal.  The  bstckground  coimt  rate 
(10)  is  determined  by 


9 


LbIrLbjLx 

he 


(2.16) 


where. 
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Lq  =  background  radiance  {Wjm^  sr) 

=  detecting  optic  axea  (m*) 

=  detector  solid  aiigle  (sr) 

1)  =  quantum  efficiency 

h  =  Planck’s  constant  { J  s) 

c  =  speed  of  light,  vax:uum  (m/s) 

A  =  average  wavelength  of  sensor  (m) 

For  a  sensor  viewing  deep  space,  the  radiance  of  the  celestial  background  \s  Lq  = 
2  X  10“®  Wjrn?  sr  (3:54).  For  a  sensor  with  a  1.6  m  diameter  collecting  optic  and  a 
CCD  detector  with  pixel  size  20  /im  (7:8),  5^  =  2.01  and  F^  w  4  x  10~^^  sr. 
Using  T)  =  .6  and  an  average  wavelength  of  A  =  600  nm,  the  count  rate  is  ia/q  = 
2,912  or  roughly  10^.  Target  count  rates  are  greater  them  background  noise  by  at 
least  an  order  of  magnitude.  (For  the  case  studies  developed  in  the  next  chapter, 
target  count  rates  range  from  10^  to  10®.)  Since  the  target  signal  dominates  the 
background  noise,  background  noise  is  assumed  negligible  in  this  analysis. 

2-4  Detector  Noise 

Detector  noise  for  a  CCD  is  from  read  noise.  The  read  noise  for  a  CCD  detector 
ranges  from  am  =  6  e~  to  am  =  20  e~  or  higher  depending  on  the  choice  of  CCD 
and  intended  application  (8).  For  this  analysis,  the  CCD  read  noise  is  assumed  to 
be  am  =  12  e“  (7),  which  corresponds  to  a  CCD  detector  having  high  quantum 
efficiency  in  the  wavelength  region  300  nm  to  900  nm  (8).  In  terms  of  a  detector 
noise  count  rate,  the  detector  noise  is  ip/q  =  12  sec~^. 

2.5  Integration  Time 

Integration  time  is  limited  by  target-sensor  viewing  geometry,  time  in  coverage, 
target  stability,  and  the  relative  velocity  between  target  and  sensor.  The  integration 
time  is  case  dependent.  A  complicating  factor  in  hyperspectral  analysis  is  the  small 
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wavelength  bandwidth.  Integration  time  must  be  longer  to  reach  the  same  SNR  as 
a  system  with  a  large  wavelength  bandwidth.  What  is  a  reasonable  integration  time 
for  hyperspectral  analysis?  The  approach  in  this  analysis  uses  a  long  integration 
time  of  =  8  s  to  determine  if  hyperspectral  anzdysis  is  feasible  for  long,  best  case 
integration  times.  If  hyperspectral  analysis  is  feasible  for  long  integration  times,  the 
shortest  possible  integration  time  is  calculated  by  setting  a  minimum  SNR  threshold 
and  solving  for  the  integration  time. 

2.6  Signal  to  Noise  Ratio 

The  general  form  of  the  SNR  equation  is  (10) 


SNR  =  .  --  ==  (2.17) 

where. 

Kg  =  average  number  of  photon  events  caused  by  desired  signal  {counts) 

Kq  =  average  number  of  photon  events  caused  by  background  {counts) 

(Trn  =  read  noise  of  detector  (e“). 

In  this  analysis,  background  is  negligible  (Section  2.3),  so  Equation  2.17  simplifies 
to 

SNR  =  ■■  .  (2.18) 

yKg+Of^ 

Equation  2.18  is  used  to  derive  SNR  equations  for  imaging  and  spectrometry  cases. 

2.6.1  SNR  Equation  for  Imaging.  For  the  imaging  case,  light  is  divided 
among  the  individual  detectors  or  pixels  in  the  detector  axray.  The  SNR  per  pixel. 
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or  SNRp^^,  is  a  measure  of  the  signal  to  noise  ratio  for  each  individual  detector. 
To  determine  if  the  sensor  is  imaging,  the  angular  subtense  of  the  diffraction  limited 
spot  is  compared  to  the  target’s  angular  subtense.  The  juigular  subtense  of  the 
diffraction  limited  spot  is 

=  1-225  (2.19) 

where, 

A  =  wavelength  (m) 

D  =  diameter  of  primary  optic  (m). 

The  target  angular  subtense  is 

^  (2-20) 


where. 


Lj  =  length  of  target  (m) 
r  =  range  to  target  (m). 

If  <  A4>sati  system  can  image  the  object.  This  is  a  rudimentary  defi¬ 

nition  of  imaging,  and  does  not  consider  resolution  criteria.  The  SNRp^^  equation, 
from  2.18  with  Kg  =  Kp^^,  is 


where. 


SNR 


K, 


■pix 


ptx 


y  ^pix  +  ^ 


(2.21) 


^pix  ~  average  number  of  photon  events  per  pixel  exposed 
<Trn  =  read  noise  of  detector  (e“). 
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To  express  K~~  in  terms  of  the  number  of  detectors  exposed  by  a  target  (10), 


Krp 

K  ■  -  ^ 

^mx  —  \T  .  2 

^^ptx 


where, 


Kj'  =  average  number  of  photon  events  caused  by  target  (counts) 
^pix  ~  ^lumber  of  pixels  exposed  by  target  in  one  dimension. 

The  average  number  of  photon  events  caused  by  the  target  (10)  is 


(2.22) 


=  iU 


(2.23) 


where, 

^  =  count  rate  from  the  detector  (1/s) 

=  integration  time  (s). 

The  number  of  pixels  exposed  by  a  target  in  one  dimension  depends  on  the  target’s 
angular  subtense  and  the  detector  angular  subtense  (10) 


(2.24) 


where, 

A<f>gaf  =  target  angular  subtense  (rad) 
das  =  detector  angular  subtense  (rad). 

The  detector  area  subtense  expressed  in  terms  of  a  detector  size  and  focal  length  is 

(10) 


(2.25) 


where. 
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d  =  detector  size  (m) 

/  =  focal  length  (m). 

To  avoid  specifying  a  focal  length  for  the  sensor,  the  detector  area  subtense  was 
assumed  to  be  equal  to  the  diffraction  limited  spot, 


das  - 


(2.26) 


When  Equations  2.19,  2.20,  2.23,  are  substituted  into  Equation  2.21,  the  SNRp^ 
expression  becomes 


SNR. 


tL 

q 


■pix 


If  Sj'  «  Lj'^  and  «  D^,  Equation  2.27  becomes 


SNR. 


pix 


q  dSj^Sji 


lir.  (1.22  rA)2 


+  <rrn^ 


(2.27) 


(2.28) 


When  the  hypemormalized  parameter  (Equation  2.12)  and  the  count  rate  equation 
(Exjuation  2.15)  are  substituted  into  the  SNRp^^  equation,  the  parameters  5^,  Sj, 
and  cancel.  When  an  object  exposes  more  than  one  detector,  imaging  is  possible, 
and  the  area  of  the  detecting  optic  5 target  area  Sj’,  and  range  r,  drop  out  of  the 
imaging  SNR  calculation.  To  make  this  approximation,  the  range,  target  area,  and 
area  of  detecting  optic  cannot  vary,  which  negates  the  use  of  the  hypemormalized 
parameter.  By  assuming  das  =  the  area  of  the  detecting  optic  is  hidden 

with  the  equation  as  part  of  the  number  of  pixels  exposed  by  the  target  Np^^.. 
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Cancelling  S Sy,  and  r^,  the  SNRp^j.  equation  becomes 


SNRj^^  « 


P  pt^  1.22^  cos(^^ )  cos(ff2)  '’’o 


yJP  p  1.22^  co8(di )  008(^2)  ‘^oV  +  (^rn  ‘ 
The  constants  and  angles  can  be  grouped  into  a  single  parameter 


(2.29) 


z  =  1.22^  cos(^j)  cos(^2)^o>7 


(2.30) 


The  parameter  2  is  a  maximum  when  the  encounter  angles  equal  zero,  and  quantum 
efficiency  is  a  maximum.  ($1  =$2=0  deg,  Tjmax  =  -7).  The  maximum  value  for 
z  is  1.02  «  1.  When  the  parameter  2  varies  between  0  zmd  1,  the  model  simulates 
encounter  angles  going  from  90  deg  to  0  deg.  The  final  expression  v  SNRpij.  is 


SNR. 


■ptx 


LeldE. 


\JP  p  A*  2  +  0‘rn  ^ 


(2.31) 


2.6.2  SNR  Equation  for  Spectrometry.  In  contrast  to  the  imaging  case,  for 
spectrometry  light  is  incident  on  a  single  detector.  The  single  pixel  SNR  is  given  by 
Equation  2.18  with  Kg  =  K'p 

SNR  =  -  (2.32) 

V^+cTrn" 

where, 

Kj>  =  average  number  of  photon  events  per  integration  time  caused  by 
target  (counts) 

arn  =  read  noise  of  detector  (e”). 
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Using  Equation  2.23  and  expressing  the  read  noise  in  terms  of  a  count  rate  from  the 
detector,  the  single  pixel  SNR  (10)  is 


where, 


SNR 


(2.33) 


k. 

9 

in 

9 

u 


=  target  signal  coimt  rate  (1/s) 
=  detector  noise  count  rate  (l/s) 
=  integration  time  (s) 


Three  tools  are  in  place  for  the  SNR  evaluation:  the  target  signal  count  rate,  and 
the  SNR  expressions  for  imaging  and  spectrometry.  The  target  signal  count  rate  was 
developed  from  basic  radiometry.  Other  elements  required  for  the  SNR  evaluation  are 
the  background  noise,  detector  noise,  and  integration  time.  Background  noise  is  small 
compared  to  the  target  signal,  and  is  assumed  negligible.  Detector  noise  is  modeled 
on  a  CCD  detector,  with  orn  =  12e~.  Integration  time  is  case  dependent,  but  can  be 
bounded  by  running  models  for  a  best  case  scenario  with  a  long  integration  time,  and 
then  setting  an  SNR  threshold  and  solving  for  the  shortest  possible  integration  time 
for  a  worst  case  scenario.  When  the  signal  is  incident  on  more  than  one  detector, 
the  SNRpi^  determines  if  enough  signal  is  present  in  each  pixel  for  imaging.  The 
single  pixel  SNR  is  useful  for  spectrometry,  when  the  diffraction  limited  spot  falls 
on  one  detector. 


III.  Results 


S.  1  Introduction 

The  target  signal  count  rate,  SNR  and  SNRp^j.  were  calculated  for  four  cases. 
Two  cases  concern  hyperspectral  imaging  and  two  concern  hyperspectrometry.  This 
chapter  outlines  the  procedures  for  running  a  case  study  and  presents  the  results  for 
the  hyperspectral  imaging  and  hyperspectrometry  case  studies. 


3.2  Procedure 

To  run  a  case  study,  the  procedure  involves  specifying  a  configuration,  deter¬ 
mining  the  spectral  reflectance  of  the  target,  running  the  radiometry  model,  and 
calculating  the  minimum  integration  time. 


3.2. 1  Configuration.  The  model  configuration  determines  the  system’s 
imaging  ability  and  the  hypernormalized  parameter.  The  imaging  determination 
compares  the  target  angular  subtense  A  (Equation  2.20)  to  the  angulau-  subtense 
of  the  diffraction  limited  spot,  (Equation  2.19).  The  A^^j^  calculation 

uses  a  m20cimum  wavelength  of  900  nm  to  give  the  largest  diffraction  limited  spot 
angular  subtense  for  the  wavelength  region  of  300-900  nm.  If  the  system  can  image 
the  target,  the  number  of  pixels,  can  be  determined  (Equation  2.24)  since 

^  ^min  = 


^^min 


The  hypernormalized  parameter  calculated  in  the  configuration  is  a  reference  point 
used  in  the  case  studies.  The  reference  hypernormalized  parameter 
maximum  fl  for  normal  incidence  of  the  target-sensor  angle  (^j)  and  the  source- 
target  angle  (Oq)-  ^ref  determines  the  bound  for  the  hypernormalized  parameter  in 
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the  radiometry  model,  and  is  the  basis  of  comparison  between  two  models.  Appendix 
A  (page  A-3)  contains  the  Mathcad  Configuration  code. 

3.2.2  Spectral  Reflectance  Interpolation.  The  spectral  reflectance  values  for 
a  given  material  are  compiled  in  the  Matter. dot  Surface  Properties  Data  Compilation 
(14).  This  compilation  contains  spectral  reflectance  values  for  materials  commonly 
used  on  spacecraft.  The  spectral  reflectance  data  is  given  in  tabular  (spectral  re¬ 
flectances  for  an  increment  of  10  nm)  and  graphical  form.  Figure  3.1  shows  an 
example  spectral  relectance  plot  for  solar  cells,  Surface  Material  Code  23  (SMC-23). 


0.0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6  1.8  2.0 

'  Wavelength  (microns) 


Figure  3.1  Spectral  Reflectance  (%),  SMC-23:  Solar  Cell,  Silicon 

Appendix  G  contains  the  spectral  reflectance  graphs  for  the  five  materials  used 
in  the  case  studies.  The  wavelength  increment  in  the  compilation  is  A  A  =  10  nm. 
To  use  the  A  A  =  20  nm  and  A  A  :=  5  nm  models,  a  Mathcad  interpolation  routine 
converts  the  10  nm  increment  into  20  nm  and  5  nm  increments.  Interpolating  to 
5  nm  assumes  that  the  spectral  reflectance  values  in  the  Matter.dat  Surface  Properties 
Data  Compilation  are  continuous.  (Appendix  A,  pages  A-4  through  A-6.) 


3-2 


S.8.S  Radiometry  Model.  The  raxliometry  model  calculates  the  count  rate, 
SNR  and  SNRpij..  Inputs  to  the  model  include: 

•  wavelength  lower  and  upper  bound 

•  wavelength  bandwidth 

•  spectral  reflect2uice  data  file 

•  range  of  the  hypernormalized  parameter 

•  integration  time 

•  detector  noise 

•  fractional  material  areas 

The  radiometry  model  outputs  data  to  graphs. 

3.8.3. 1  Wavelength  Lower  and  Upper  Bound.  From  Section  2. 2. 2.2, 
the  wavelength  region  is  300-900  nm.  To  match  the  spectral  reflectance  values  to  the 
center  of  the  wavelength  bandwidth,  =  300  — (A  A/2)  and  A^^j  =  900  + (A  A/2). 
The  wavelength  lower  and  upper  bound,  the  wavelength  bandwidth,  eind  a  counter 
(i)  are  used  to  determine  the  limits  of  integration  for  the  photon  flux  calculation. 

3. 8. 3. 8  Wavelength  Bandwidth.  A  A  =  20  nm,  10  nm,  or  5  nm. 

3. 8.3. 3  Spectral  Reflectance  Data  File.  The  spectral  reflectance  is 
read  from  a  data  file  specific  to  material  and  wavelength.  The  term  P23  means 
spectral  reflectance  for  Surface  Material  Code  23  (SMC-23). 

3. 8. 3. 4  Hypernormalized  Parameter.  The  hypernormalized  parame¬ 
ter  adlows  the  modeling  of  a  wide  range  of  cases  in  one  radiometry  calculation.  Section 
3. 3. 1.5  illustrates  the  application  of  for  Case  (1).  To  eiccomplish  the  modeling, 
the  bound  of  the  hypernorm£dized  parameter  Q  is  varied  by  an  order  of  magnitude 
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above  and  below  exponent.  For  example,  if  =  10  the  reinge  for  fi 

is  10~12 

5.2.3.5  Integration  Time.  From  Section  2.5,  the  integration  time  is 
t^  =  8  s.  This  models  a  best  case  scenario.  Section  3.2.4  discusses  the  best  and 
worst  case  integration  time  approach  for  this  analysis. 

5.2.5.6  Detector  Noise.  From  Section  2.4,  the  detector  noise  is  Z£)  = 
12 9  scc~^  or  ffrn  =  12. 

5.2.5.7  Fractional  Material  Areas.  When  a  case  study  involves  an 
object  made  of  different  materials,  the  fractional  material  areas  are  used  to  weight 
the  target’s  count  rate.  For  example,  an  object  made  of  70  percent  solar  array 
(SMC- 23)  and  30  percent  aluminum  (SMC-02)  would  have  a  count  rate  expression 
weighted  as 


'T{total)  =  •’^0i2’(23)  +-30*r(02) 

3.2.4  Example  Input  and  Output.  With  all  the  input  parameters  defined, 
the  radiometry  model  executes  and  presents  the  results  in  three  graphs.  An  example 
of  the  input  required  and  resulting  output  follows.  The  excimple  is  for  Case  (1), 
AA  =  10  nm. 

Input  parameters: 

•  wavelength  lower  and  upper  bound:  A^j  =  295  nm  and  A^j  =  905  nm 

•  wavelength  bandwidth:  A  A  =  10  nm 

•  spectral  reflectance  data  file:  smc02xx.dat  {smc02=Surface  Material  Code  02; 
XX  for  20  nm) 

•  range  of  the  hypemormalized  parameter:  fl  =  10~^^  . . .  10~^^ 
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•  integration  time:  =  Ss 

•  detector  noise:  ijj  =  \2qstc~^  or  arn  =  12 

•  fractional  material  areas:  None  (object  is  maule  entirely  of  SMC-02) 

Output:  The  output  from  the  radiometry  model  is  three  graphs;  Figure  3.2,  Fig¬ 
ure  3.3,  and  Figure  3.4^.  Figure  3.2  shows  the  count  rate  {ij'lq)  vs  fl  for  Case  (1), 
AA  =  10 nm.  The  graph  is  a  log-log  plot,  with  the  target  signal  count  rate,  ij'lq, 
graphed  on  the  vertical  axis,  and  0  graphed  on  the  horizontal  axis.  The  vertical 
dashed  line  is  the  reference  hypernormalized  p^u^ameter,  computed  for  the 

target.  The  legend  for  converting  the  counter(t)  to  the  mean  wavelength  of  the  in¬ 
terval  X  (Appendix  A,  page  A-8)  is  given  in  Table  3.1.  This  legend  applies  to  all  three 


Table  3.1  Legend,  Case  (1),  10  nm 

graphs.  Figure  3.3  shows  the  SNR  vs  fl  for  Case  (1),  A  A  =  10  nm.  The  graph  uses 
a  log  plot  for  the  U  axis  and  normal  plot  for  the  SNR  on  the  vertical  axis.  The  verti¬ 
cal  dashed  line  is  the  reference  hypernormalized  parameter,  ^ref'  computed  for  the 
target.  Figure  3.4  shows  the  SNRp^j.  vs  z  for  Case  (1),  AA  =  10 nm.  SNRp^j.  is 
plotted  on  the  vertical  axis,  and  z,  the  cingles  parameter,  on  the  horizontal  axis.  The 
parameter  z  varies  from  0-1,  simulating  encounter  angles  9i  and  $2  going  from  90 
deg-0  deg.  All  graphs  in  this  chapter  and  the  appendices  follow  these  three  formats. 

3.2.5  Integration  Time.  The  integration  time  for  this  analysis  considers 
a  best  case  and  worst  case  situation.  The  best  case  is  for  a  stable  target  slowly 

^Figures  used  in  this  chapter  are  duplicated  in  the  appendix. 


3-5 


9.494*10 


Figure  3.3  Case  (1),  SNR  vs  ft{radm^),  10  nm 
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Figure  3.4  Case  (1),  SNR—j.  vs  z,  10  nm 
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crossing  the  field  of  view  of  the  sensor,  with  little  aspect  change.  For  this  situation, 
the  integration  time  is  set  at  =  8  s,  the  value  used  in  the  radiometry  models.  The 
worst  case  is  when  the  target-sensor  relative  velocity  is  high,  or  the  target  is  tumbling, 
causing  radical  aspect  changes.  In  this  case,  the  integration  time  required  to  achieve 
a  minimum  SNR  or  SNRp^j.  value  was  obtained  by  solving  for  the  integration  time 
in  Equation  2.31  and  Equation  2.33.  When  imaging  is  not  possible,  the  equation  for 
the  integration  time  is 


u  = 


(3.1) 


When  imaging  is  possible,  the  equation  for  the  integration  time  is 


‘d  = 


SNR^^^(zH)±  +4(2g)^  SNRj^^Urn' 


2{zHy 


(3.2) 


where, 

H 

z  =  1.22^  cos(^l)  cos{02)  To  rj 

For  this  analysis,  the  minimum  SNRp^^  and  SNR  thresholds  are  three. 

S.3  Hyperspectral  Imaging 

The  hyperspectral  imaging  case  studies  are  for  a  rocket  body  in  low  earth 
orbit.  The  case  studies  are  identical,  except  for  Case  (1)  having  an  aluminum  rocket 
body  and  Case  (2)  having  the  rocket  body  covered  with  white  paint  (9).  After 
the  individual  case  results  are  discussed,  the  two  cases  are  compared  with  each 
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other.  Three  key  elements  in  this  comparison  are  signal  strength,  line  order,  and 
line  separation. 

S.3.1  Case  One. 

3. 3. 1.1  Configuration.  Target:  The  modeled  rocket  body  is  6  m  long 
and  2  m  in  diameter.  The  range  to  the  target  is  zissumed  to  be  800  km. 

Sensor:  The  sensor  is  a  1.6  m  diameter  telescope.  The  mean  wavelengths  of  interest 
are  450  nm,  550  nm,  650  nm,  750  nm  and  850  nm^. 

Imaging  determination:  The  maximum  target  angular  subtense  is 

A<l>gat  =  7.500  X  10“^  rad 

and  the  angular  subtense  of  the  diffraction  limited  spot  is 

A<j>min  —  6.863  X  10~^  rad. 

A4>min  ^^sah  so  imaging  is  possible.  The  number  of  pixels  exposed  by  the 
target  is  «  11, 

HypernormaJized  parameter:  The  reference  hypernormalized  parameter  for  normal 
incidence  =^2  =  0  deg)  is  (Equation  2,12) 

=  2.217  X  10~^^  radm^ 

The  range  of  the  hypernormalized  parameter  in  the  radiometry  model  is 
n  =  10“12...io-10radm2.. 

3. 3. 1.2  Spectral  Reflectance.  The  rocket  body  is  made  of  aluminum, 
SMC-02  (Figure  F.l). 

^For  AA  =  20  nm,  an  increment  of  20  nm  beginning  at  300  nm  will  only  give  wavelength 
intervals  of  interest  at  440  nm,  540  nm,  640  nm,  740  nm,  and  840  nm. 
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3.3. 1.3  Radiometry  results.  Table  3.2,  Summary  of  Results,  Case 
(1),  shows  enough  signal  is  present  in  each  of  the  wavelength  bandwidths  to  perform 
hyperspectral  imaging.  The  count  rate  and  SNR  values  in  Table  3.2  correspond  to 
the  points  where  i^itersects  the  plotted  lines  in  Figures  B.1-B.9,  giving  MIN 
and  MAX  values  as  a  function  of  mean  wavelength.  The  SNRj^^  values  correspond 
to  value  z  —  1  (0^  =  ^2  =  0  deg).  The  lowest  count  rate  is  9.5  x  10^,  which  is 
roughly  five  orders  of  magnitude  above  the  background  coimt  rate  determined  in 
Section  2.3.  The  lowest  values  of  SNR  and  SNRp^^.  are  well  above  the  minimum 
threshold  of  three.  Appendix  B  contains  a  complete  set  of  graphs  for  Case  (1), 


CASE  1 
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gth  Bandwidth  (A  A) 
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Table  3.2  Summary  of  Results,  Case  (1) 


Figure  B.l  to  Figure  B.9.  The  set  consists  of  a  legend  (Table  B.l)  converting  i  to  A 
and  graphs  of  count  rate,  SNR,  and  SNRp^j.  for  each  of  the  20  nm,  10  nm,  and  5 
nm  wavelength  bandwidths.  The  10  nm  and  5  nm  curves  are  very  similar  due  to  the 
interpolation  of  the  spectral  reflectance  values  (compare  Figure  B.5  to  Figure  B.8). 

The  line  order  and  spacing  show  the  target  signal’s  dependence  on  wavelength 
and  bandwidth.  The  order  of  the  lines  is  determined  by  starting  at  the  lowest 
plotted  line  (bottom)  corresponding  to  the  mean  wavelength  with  the  weakest  target 
signal,  and  listing  successive  lines  until  reaching  the  line  corresponding  to  the  mean 
wavelength  with  the  largest  target  signal  (top).  (Signal  strength,  line  order,  and 
line  spacing  are  used  as  the  basis  of  comparison  between  targets.)  Line  order  is 
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ORDER 

20  nm 

10  nm 

5  nm 

Top  5 

640 

650 

650 

4 

540 

550 

550 

3 

740 

750 

750 

2 

840 

850 

850 

Bottom  1 

440 

450 

450 

Missing  honz<^tal  lines  indicate 
overlapping  wavelength  intervals 


Table  3.3  Line  Order,  Case  (1) 

not  important  for  S  since  SNRj^^,  is  not  dependent  on  target  area,  area  of 

detecting  optic,  and  range  (see  Section  2.6.1).  The  line  order  for  Case  (1),  Table  3.3, 
does  not  change  for  the  three  wavelength  bandwidths.  Two  lines  overlap  for  the 
20  nm  bandwidth,  the  540  nm-740  nm  pair  (Figme  B.2).  Overlapping  lines  do  not 
have  a  line  separating  the  two  wavelength  intervals.  This  convention  is  used  in  all 
line  order  tables.  Possible  reasons  for  the  difference  between  the  20  nm  and  other 
cases  are  the  wider  wavelength  interval  and  coarser  spectral  reflectance  values  from 
the  interpolation  of  the  spectral  reflectance  data  file. 

3. 3.1.4  Integration  time.  Since  imaging  is  possible,  the  minimum 
integration  times  are  calculated  using  Equation  3.2.  The  minimum  integration  times 
for  SNRpi^  =  3  for  each  wavelength  bandwidth  are  given  in  Table  3.4.  As  the 
bandwidth  decreeises,  the  integration  time  increases.  An  integration  time  in  mi¬ 
croseconds  is  very  small;  however,  encoimter  mechanics  such  as  orbital  velocity  and 
target  stability  may  necessitate  a  very  small  integration  time. 


AA 

td 

20  nm 

12.8  fisec 

10  nm 

23.5  fisec 

5  nm 

47.1  /isec 

Table  3.4  Minimum  integration  time.  Case  (1) 


S.S.1.5  Using  the  hypemormalized  parameter  fi.  Why  use  a  hy- 
pernormalized  parameter?  The  advantage  is  one  graph  indicates  a  wide  range  of 
target-sensor  encounters.  Three  examples  illustrate  the  use  of  the  hypemormalized 
parameter  fi:  varying  encoimter  angles;  varying  target  range;  and  varying  the  diam¬ 
eter  of  the  sensor’s  collecting  optic.  A  requirement  for  these  variations  is  to  remain 
within  the  bounds  of  fl  for  Case  (1),  10“^^  <  fl  <  The  illustrations  will 

use  Figure  3.3,  Case  1,  SNR,  10  nm  to  show  how  varying  (I  applies  to  the  graphs. 
Since  the  illustrations  use  the  single  pixel  SNR,  the  imaging  determination  will  not 
be  recalculated  for  varying  range  and  collecting  optic  diameter. 

Basic  Case  (1):  For  comparison  purposes.  Case  (l)’s  ^ref  =  $2  =  0  deg  is 

(Equation  2.12) 


^ref  ~  2.217x10  ^^radm^. 

Single  pixel  SNR  values  from  Table  3.2  are  SNRmax  =  4.4  x  10^  and  SNRmin  = 
3.9  X  10^. 

Varying  encounter  angles:  As  the  source-target  angle  and  target-sensor  angle  ap¬ 
proach  90  deg,  the  SNR  will  decrease.  The  hypemormalized  parameter  for  $1  = 
75  deg  and  $2  =  75  deg  is 


n  =  1.485  x  10  ^^radrn^. 

Single  pixel  SNR  values  from  Figure  3.3  are  SNRmax  =  1.0  x  10^  and  SNRjam  = 
0.9  X  10^. 

Varying  range:  As  the  range  increases,  the  SNR  will  decrease.  The  hypemormalized 
parameter  for  r  =  3000  km  is 

n  =  1.576  X  10“^2^q^^2_ 
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Single  pixel  SNR  values  from  Fig-'re  3.3  are  SNRmax  =  1.0  x  10^  emd  SNRmjn  = 
0.9  X  lO'^. 

Varying  sensor  collecting  optic  diameter:  As  the  diameter  decreases,  the  SNR  will 
decrease.  The  hypemormalized  parameter  for  D  =  .2  m  and  r  =  800  km  is  Q  = 
3.464  X  10~^^  rod m^,  which  is  outside  the  bound  of  fi.  For  a  shorter  range,  r  = 
440  km,  the  hypemormalized  parameter  is 

Q  =  1.385  X  10~^^  rad  m^. 

Single  pixel  SNR  values  from  Figure  3.3  are  SNRmax  =  1.0  x  10^  and  SNRmin  = 
0.9  X  10'*. 

Each  of  the  three  different  illustrations  dropped  the  hypemormalized  parameter  to 
about  the  same  value,  fl  «  1.5  x  lO”"*^.  The  hypemormalized  parameter  was 
computed  using  the  configuration  code  (Appendix  A,  page  A-3),  a  quick  calculation 
compared  to  the  computationly  intensive  radiometry  model.  Without  mnning  the 
radiometry  model  again,  the  hypemormalized  parameter  modeled  a  variety  of  target¬ 
sensor  encounters. 

3.3.2  Case  Two. 

3.3.2. 1  Configuration.  The  configuration  for  Case  (2)  is  identical  to 

Case  (1). 

3. 3. 2. 2  Spectral  Reflectance.  The  modeled  rocket  body  is  painted 
white,  SMC-20,  instead  of  having  a  bare  aluminum  surface  (Figure  F.2). 

3. 3. 2. 3  Radiometry  results.  Table  3.5,  Summary  of  Results,  Case(2), 
shows  enough  signal  is  present  in  each  of  the  wavelength  bandwidths  to  perform 
hyperspectral  imaging.  The  count  rate  and  SNR  values  in  Table  3.5  correspond  to 
the  points  where  intersects  the  plotted  lines  in  Figures  C.1-C.9,  giving  MIN 
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and  MAX  values  as  a  function  of  mean  wavelength.  The  SNRp^j.  values  correspond 
to  value  z  =  1  {$i  =  $2  =  0  deg).  The  lowest  count  rate  is  1.3  x  10^,  which  is 
five  orders  of  magnitude  above  the  background  count  rate  determined  in  Section  2.3. 
The  lowest  values  of  SNR  and  SNRp^^  are  well  above  the  minimum  threshold  of 
three.  Appendix  C  contains  a  complete  set  of  graphs  for  Case  (2),  Figure  C.l  to 


CASE  2 

Wavelen 

gth  Bandwidth  (A  A) 

20  nm 

10  nm 

5  nm 

Count  rate 

MAX 

msuiji 

MIN 

QQQH 

SNR 

MAX 

EQQ1E9 

QQQOm 

MIN 

MAX 

MIN 

ItfkilHikl 

iiKiiiniiki 

Table  3.5  Summary  of  Results,  Case  (2) 


Figure  C.9.  The  set  consists  of  a  legend  (Table  C.l)  converting  i  to  A  and  graphs  of 
count  rate,  SNR,  and  SNRp^^  for  each  of  the  20  nm,  10  nm  and  5  nm  wavelength 
bandwidths.  The  10  nm  and  5  nm  curves  are  very  similar  due  to  the  interpolation  of 
the  spectral  reflectance  values  (compare  Figure  C.5  to  Figure  C.8).  The  line  order 
for  Case  (2),  Table  3.6,  does  not  change  for  the  three  wavelength  bandwidths.  Two 
lines  overlap  for  the  10  nm  and  5  nm  bandwidth,  the  650  nm-550  nm  pair.  Possible 
reasons  for  the  slight  difference  between  the  20  nm  and  10  nm  plots  are  the  wider 
wavelength  interval  and  coarser  spectral  reflectance  values  from  the  interpolation  of 
the  spectral  reflectance  data  file. 

S.S.2.4  Integration  time.  The  minimum  integration  times  for  each 
wavelength  bandwidth  for  an  SNRp^^  =  3  are  given  in  Table  3.7.  As  the  bandwidth 
decreases,  the  integration  time  increases.  An  integration  time  in  microseconds  is  very 
small;  however,  encounter  mechanics  such  as  orbitad  velocity  and  target  stability  may 
necessitate  a  very  small  integration  time. 
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ORDER 

20  nm 

10  nm 

5  nm 

Top  5 

640 

650 

4 

540 

550 

3 

740 

750 

2 

840 

850 

Bottom  1 

440 

450 

Missing  horizontal  lines  indicate 
overlapping  wavelength  intervals 


Table  3.6  Line  Order,  Case  (2) 


AA 

td 

20  nm 

9.6  /isec 

10  nm 

17.7  psec 

5  nm 

35.4  psec 

Table  3.7  Minimum  integration  time,  Case  (2)  . 

3.3.3  Comparison,  Case  One  versus  Case  Two.  For  this  analysis,  the  three 
areas  of  comparison  are  the  signal  strength  from  the  target,  the  line  order,  and  line 
spacing.  This  comparison  involves  Case  (1)  and  Case  (2)  single  pixel  SNR  graphs 
with  AA=  10  nm.  Legends  for  the  Case  (1)  and  Case  (2)  graphs  are  the  same,  and 
given  in  Table  3.1.  Case  (1)  SNR  10  nm  was  used  earlier  in  Section  3.2.3. 8,  (see 
Figure  3.3).  The  Case  (2)  counterpart.  Case  (2)  SNR  10  nm  is  shown  in  Figure  3.5. 


ORDER 

20 

nm 

10  nm 

5  nm 

Case  (1) 

Case  (2) 

Ceae  (2) 

Case  (2) 

Top  5 

640 

640 

650 

650 

MB 

4 

540 

550 

kb 

550 

KB 

3 

740 

750 

750 

750 

750 

2 

840 

840 

850 

850 

850 

Bottom  1 

440 

*  .  1  4«  ' 

440 

450 

450 

- - 7T 

450 

450 

Table  3.8  Line  Order,  Case  (1)  and  Case  (2)  Combined 
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To  aid  the  comparison,  the  line  orders  of  Case  (1)  and  Case  (2)  are  combined 
in  Table  3.8.  The  line  spacing  comparison  is  made  using  Table  3.8,  Figure  3.3,  and 
Figure  3.5. 

Area  Comparison 

Signal  Strength  Case  (2)  has  larger  SNR  values  than  Case  (1). 

Line  Order  Line  order  is  the  same  for  the  two  cases. 

Line  Spacing  Line  spacing  is  different  for  all  bandwidths. 

Space  Object  Identification  information  may  be  stored  in  the  different  target  signal 
strengths  in  different  spectral  bands.  The  differences  between  Case  (1)  and  Case  (2) 
are  from  the  different  spectral  reflectances  of  aluminum  and  white  paint. 

3.4  Hyperspectrometry 

The  two  cases  for  hyperspectrometry  involve  geosynchronous  communications 
satellites.  The  satellites  are  cylindrical  with  an  antenna  on  top,  but  have  different 
sizes  (9).  The  results  for  each  case  are  discussed  individually,  and  then  compared  to 
each  other. 

3.4- 1  Case  Three. 

3. 4-1.1  Configuration.  Target:  A  cylindrical  communications  satel¬ 
lite  is  7.6  m  long  and  3.64  m  in  diameter.  The  satellite  has  an  antenna  2.4  m  in 
diameter  on  top  of  the  cylinder  (1:12).  The  range  to  the  taxget  is  36,000  km,  the 
altitude  of  a  geosynchronous  satellite. 

Sensor:  The  sensor  is  a  1.6  m  diameter  telescope.  Three  regions  of  the  electromag¬ 
netic  spectrum  are  of  interest:  300  nm-400  nm,  Ultraviolet  (UV);  380  nm-540  nm. 
Blue  (B);  and  480  nm-650  nm,  Visible  (V).  Within  each  region,  the  mean  wave¬ 
lengths  of  interest  are  in  increments  of  20  nm.  (9) 
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Imaging  determination:  The  maximum  target  angular  subtense  is 

=  2.778  X  rad 

and  the  angular  subtense  of  the  diffraction  limited  spot  is 

^^min  ~  6.863  x  10~^  rad 
^  ^^sah  ^  imaging  is  not  possible. 

Hypernormalized  parameter:  The  reference  hypernormalized  parameter  for  normal 
incidence  =  ^2  =  0  deg)  is 

=  2.936  X  10“^^  rad 

The  range  of  the  hypernormalized  parameter  in  the  radiometry  model  is 
(I  =  10”^^  . . .  10~^^  radm^. 

3. 4-1. 2  Spectral  Reflectance.  The  satellite  was  modelled  with  three 
materials:  solax  cells,  SMC-23  (Figure  F.4);  a  mirrored  thermal  control  band  (as¬ 
sumed  to  be  silver),  SMC-18  (Figure  F.3);  and  aluminized  mylar,  SMC-30  (Fig¬ 
ure  F.5).  The  thermal  control  band  is  approximately  1/7  the  length  of  the  cylinder. 
The  fractional  material  areas  are  73.7  percent  solax  cells,  12.3  percent  silver,  and 
14.0  percent  aluminized  mylax,  and  yield  a  count  rate  expression  of 

*T{total)  =  •'^3’^*r(23)  -|-.140ij^30). 

3.4. 1.3  Radiometry  results.  Table  3.9,  Summary  of  Results,  Case(3), 
shows  enough  signal  is  present  in  each  of  the  wavelength  bandwidths  to  perform 
hyperspectrometry.  The  count  rate  and  SNR  values  in  Table  3.9  correspond  to  the 
points  where  intersects  the  plotted  lines  in  Figures  D.1-D.18,  giving  MIN  and 


3-19 


MAX  values  as  a  function  of  mean  wavelength.  The  lowest  count  rate  is  1.9  x  10^, 
which  is  an  order  of  magnitude  above  the  background  count  rate  determined  in 
Section  2.3.  The  lowest  value  of  the  SNR  is  well  above  the  minimum  threshold  of 
three. 


CASE  3 

Wavelength  Bandwidth  (A  A) 

20  nm 

10  nm 

5  nm 

Count  rate 

UV 

MAX 

2.5  X  10^ 

1.3  X  10^ 

6.4  X  10^ 

MIN 

7  5  X  10^ 

3.7  X  lO'^ 

1.9  X  10'^ 

B 

MAX 

2.7  X  10^ 

1.4  X  10^ 

5.4  X  lO'* 

MIN 

7.8  X  10^ 

8.5  X  10^ 

3.2  X  10^ 

V 

MAX 

2.9  X  10^ 

1.5  X  10^ 

7.0  X  10^ 

MIN 

1.1  X  10^ 

9.5  X  lO'^ 

4.9  X  lO'^ 

SNR 

UV 

MAX 

1.4  X  10^ 

0.99  X  10^ 

0.73  X  10^ 

MIN 

0.75  X  10^ 

0.54  X  10^ 

0.38  X  10^ 

B 

MAX 

1.5  X  10^ 

1.0  X  10^ 

0.71  X  10^ 

MIN 

1.1  X  10^ 

0.83  X  10^ 

0.58  X  10^ 

V 

MAX 

1.5  xlO^ 

1.0  X  10^ 

0.75  X  10^ 

MIN 

1.2  X  10^ 

0.88  X  10^ 

0.62  X  10^ 

Table  3.9  Summary  of  Results,  Case  (3) 


Appendix  D  contains  a  complete  set  of  graphs  for  Case  (3),  Figure  D.l  to 
Figure  D.18.  The  set  consists  of  a  legend  (Table  D.l)  converting  i  to  A  and  graphs  of 
count  rate  and  SNR  for  each  of  the  20  nm,  10  nm  and  5  nm  wavelength  bandwidths. 
The  10  nm  and  5  nm  curves  are  very  similar  due  to  the  interpolation  of  the  spectral 
reflectance  vzJues  (compare  Figure  D.IO  to  Figure  D.16).  The  line  order  for  Case 
(3)  is  separated  into  Ultraviolet  (UV),  Blue  (B)  and  Visible(V)  regions. 

Table  3.10  shows  the  line  order  for  the  UV  region.  The  20  nm  bandwidth 
has  a  different  order  than  the  10  nm  and  5  nm.  Possible  reasons  for  this  difference 
are  the  wider  wavelength  interval  and  coarser  spectral  reflectance  values  from  the 
interpolation  of  the  spectral  reflectance  data  file.  Two  lines  overlap  for  the  10  nm 
and  5  nm  beindwidth,  the  380  nm-360  nm  paur. 
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ULTRAVIOLET  (UV 

20  nm 

10  nm 

5  nm 

360 

380 

5 

380 

360 

360 

4 

340 

400 

3 

400 

340 

340 

2 

320 

320 

320 

Bottom  1 

m — r—  1 

300 

300 

300 

Missing  horizontal  lines  indicate 
overlapping  wavelength  intervals 


Table  3.10  Line  Order,  Case  (3),  UV 


BLUE  (B) 


20  nm 

10  nm 

5  nm 

Top  9 

540 

540 

540 

8 

520 

380 

380 

7 

|||Kg« 

520 

520 

6 

400 

400 

5 

480 

500 

500 

4 

400 

480 

3 

460 

420 

420 

2 

420 

r  460 

460 

Bottom  1 

440 

440 

440 

Missing  horizontal  lines  indicate 
overlapping  wavelength  intervals 


Table  3.11  Line  Order,  Case  (3),  B 


Table  3.11  shows  the  line  order  for  the  B  region.  The  20  nm  bandwidth  has 
a  dilFerent  order  than  the  10  nm  and  5  nm.  Possible  reasons  for  this  difference 
are  the  wider  wavelength  interval  and  coarser  spectral  reflectance  values  from  the 
interpolation  of  the  spectral  reflectance  data  file.  Two  line  pairs  overlap:  the  380 
nm-540  nm  pair  for  the  10  nm  and  5  nm  bandwidths,  and  the  380  nm-500  nm  pair 
for  the  20  nm  bandwidth. 


VISIBLE  (V) 

ORDER 

20  nm 

10  nm 

5  nm 

Top  9 

640 

640 

640 

8 

620 

620 

620 

7 

600 

600 

600 

6 

580 

580 

580 

5 

560 

560 

560 

4 

540 

540 

540 

3 

520 

520 

520 

2 

500 

500 

500 

Bottom  1 

480 

480 

480 

Missing  horizontal  lines  indicate 
overlapping  wavelength  intervals 

Table  3.12  Line  Order,  Case  (3),  V 


Table  3.12  shows  the  line  order  for  the  V  region.  The  line  order  is  the  same  for 
all  three  bandwidths.  Three  lines  overlap,  those  for  600  nm,  620  nm,  and  640  nm. 

3.4-  1-4  Integration  time.  Since  imaging  is  not  possible.  Equation  3.1 
was  used  to  calculate  the  minimum  integration  time.  The  minimum  integration 
times  for  an  SNR  =  3  for  each  wavelength  bandwidth  are  given  in  Table  3.13.  As 
the  bandwidth  decreases,  the  integration  time  increases. 

3.4-2  Case  Four. 

3. 4- 2.1  Configuration.  Target:  A  cylindrical  communications  satel¬ 
lite  is  4.77  m  long  and  2.2  m  in  diameter.  The  satellite  has  an  antenna  1.83  m  in 
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20  nm  »  125  ^sec 
10  nm  248  nsec 
5  nm  496  /isec 

Table  3.13  Minimum  integration  time,  Case  (3) 

diameter  on  top  of  the  cylinder  (1:12).  The  range  to  the  target  is  36,000  km,  the  alti¬ 
tude  of  a  geosynchronous  satellite.  Sensor:  The  sensor  is  a  1.6  m  diameter  telescope. 
Three  regions  of  the  electromagnetic  spectrum  are  of  interest:  300  nm-400  nm.  Ul¬ 
traviolet  (UV);  380  nm-540  nm.  Blue  (B);  and  480  nm-650  nm.  Visible  (V).  Within 
each  region,  the  mean  wavelengths  of  interest  are  in  increments  of  20  nm  (9). 
Imaging  determination:  The  maximum  target  angular  subtense  is 

^^sat  ~  1*833  X  10~^  rad 

and  the  angular  subtense  of  the  diffraction  limited  spot  is 

A<f)min  ~  6.863  x  10“^  rad 

A<f>min  >  ^^sati  imaging  is  not  possible. 

HypemormaJized  parameter:  The  reference  hypernormalized  parameter  for  normal 
incidence  =  ^2  =  0  deg)  is 

^ref  ~  1*197  X  10~^^  rodm^ 

The  range  of  the  hypernormalized  parameter  in  the  radiometry  model  is 
n  =  rad  m^. 


3. 4-2. 2  Spectral  Reflectance.  The  satellite  was  modelled  with  three 
materials:  solar  cells,  SMC- 23  (Figure  F.4);  a  mirrored  thermal  control  band  (as- 
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sumed  to  be  silver),  SMC-18  (Figure  F.3);  and  aluminized  mylar,  SMC-30  (Fig¬ 
ure  F.5).  The  thermal  control  band  is  approximately  1/7  the  length  of  the  cylinder. 
The  fractional  material  areas  are  68.7  percent  solax  cells,  11.5  percent  silver,  and 
19.8  percent  aluminized  mylar,  and  yield  a  count  rate  expression  of 


^Titotal)  =  •687t2’(23)  +-11%(18)  +•^^^*>(30) 

S.4-2.S  Radiometry  results.  Table  3.14,  Summary  of  Results,  Case(4), 
shows  enough  signal  is  present  in  each  of  the  wavelength  bandwidths  to  perform  hy¬ 
perspectrometry.  The  count  rate  and  SNR  values  in  Table  3.9  correspond  to  the 
points  where  ^ref  i^^fersects  the  plotted  lines  in  Figures  D.1-D.18,  giving  MIN  and 

O 

MAX  values  as  a  function  of  mean  wavelength.  The  lowest  count  rate  is  8.0  x  10*^, 
which  is  roughly  an  order  of  magnitude  above  the  background  count  rate  determined 
in  Section  2.3.  The  lowest  value  of  the  SNR  is  well  above  the  minimum  threshold  of 
three. 


CASE  4 

Wavelen 

gth  Bandwidth  (A  A) 

20  nm 

10  nm 

5  nm 

Count  rate 

UV 

MAX 

1.0  X  10^ 

5.5  X  10^ 

2.8  X  lO'^ 

MIN 

3.2  X  10'^ 

1.7  X  10^ 

8.0  X  10^ 

B 

MAX 

1.3  X  10^ 

6.0  X  10^ 

3.0  X  lO'^ 

MIN 

7.9  X  lO'^ 

4.0  X  lO'^ 

2.0  X  10^ 

V 

MAX 

1.4  X  10^ 

6.7  X  10^ 

3.2  X  10^ 

MIN 

9.5  X  10^ 

4.7  X  10^ 

2.4  X  lO'^ 

SNR 

UV 

MAX 

0.89  X  10^ 

0.66  X  10^ 

0.47  X  10^ 

MIN 

0.52  X  10^ 

0.36  X  10^ 

0.26  X  10^ 

B 

MAX 

0.90  X  10'^ 

0.69  X  lO"^ 

0.49  X  10^ 

MIN 

0.76  X  10^ 

0.56  X  lO"^ 

0.41  X  10^ 

V 

MAX 

1.0  X  lO'^ 

0.72  X  10^ 

0.51  X  10^ 

MIN 

0.85  X  10^ 

0.60  X  10^ 

0.43  X  10^ 

Table  3.14  Summary  of  Results,  Case  (4) 
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Appendix  E  contains  a  complete  set  of  graphs  for  Case  (4),  Figure  D.l  to 
Figure  D.18.  The  set  consists  of  a  legend  (Table  E.l)  converting  t  to  A  and  graphs  of 
count  rate  and  SNR  for  each  of  the  20  nm,  10  nm,  and  5  nm  wavelength  bandwidths. 
The  10  nm  and  5  nm  curves  are  very  similar  due  to  the  interpolation  of  the  spectral 
reflectance  values  (compare  Figure  E.IO  to  Figure  E.16).  The  line  order  for  Case  (4) 
is  separated  into  Ultraviolet  (UV),  Blue  (B)  and  Visible(V)  regions. 


ULTRAVIOLET  (UV) 

ORDER 

20  nm 

10  nm 

5  nm 

TOP  6 

380 

380 

5 

380 

360 

360 

4 

400 

400 

3 

340 

340 

340 

2 

mm 

320 

320 

BOTTOM  1 

300 

300 

300 

Missing  horizontal  lines  indicate 
overlapping  wavelength  intervals 

Table  3.15  Line  Order,  Case  (4),  UV 


Table  3.15  shows  the  line  order  for  the  UV  region.  The  20  nm  bandwidth 
has  a  different  order  than  the  10  nm  and  5  nm.  Possible  reasons  for  this  difference 
are  the  wider  wavelength  interval  and  coarser  spectral  reflectance  values  from  the 
interpolation  of  the  spectral  reflectance  data  file.  One  line  pair  overlaps  for  the  20  nm 
bandwidth,  the  400  nm-340  nm  pair. 

Table  3.16  shows  the  line  order  for  the  B  region.  The  20  nm  bandwidth  has 
a  different  order  than  the  10  nm  and  5  nm.  Possible  reasons  for  this  difference 
are  the  wider  wavelength  interval  and  coarser  spectral  reflectance  values  from  the 
interpolation  of  the  spectral  reflectance  data  file.  Two  line  pairs  overlap  for  the 
20  nm  bandwidth,  the  440  nm-420  nm  pair  and  the  400  nm-460  nm  pair.  One  line 
pair  overlaps  for  the  10  nm  and  5  nm  bandwidths,  the  460  nm-420  nm  pair. 

Table  3.17  shows  the  line  order  for  the  V  region.  The  line  order  is  the  same 
for  all  three  bandwidths.  Three  lines  overlap,  the  600  nm-620  nm-640  nm. 
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issing  horizontal  lines  indicate 
overlapping  wavelength  intervals 

Table  3.16  Line  Order,  Case  (4),  B 


Missing  horizontal  lines  indicate 
overlapping  wavelength  intervals 

Table  3.17  Line  Order,  Case  (4),  V 
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S.4-2.4  Integration  time.  The  minimum  integration  times  for  an 
SNR  =  3  for  each  wavelength  bandwidth  are  given  in  Table  3.18.  As  the  bandwidth 
decreases,  the  integration  time  increases. 


AA 

u 

20  nm 

TIB  fisec 

10  nm 

553  psec 

5  nm 

1000  /usee 

Table  3.18  Minimum  integration  time,  Case  (4) 


3-4 -3  Comparison,  Case  Three  versus  Case  Four.  The  three  areas  of 
comparison  are  the  signal  strength  £rom  the  target,  the  line  order,  and  line  spacing. 
This  comparison  involves  Case  (3)  and  Case  (4)  single  pixel  SNR  graphs  with  AA= 
10  nm.  Legends  for  the  Case  (3)  and  Case  (4)  graphs  are  the  same,  and  given  in 
Table  3.19.  Comparisons  for  Case  (3)  and  Case  (4)  are  broken  into  Ultraviolet,  Blue, 
and  Visible  categories.  To  aid  the  comparison,  the  line  orders  from  Case  (3)  and 
Case  (4)  are  combined  within  each  region.  The  line  order  tables  contain  all  three 
wavelength  bandwidths. 


UV 

B 

V 

i 

A  (nm) 

{  '  A  (nm) 

i 

A  (nm) 

1 

300 

9 

380 

19 

480 

3 

320 

11 

400 

21 

500 

5 

340 

13 

420 

23 

520 

7 

360 

15 

440 

25 

540 

9 

380 

17 

460 

27 

560 

11 

400 

19 

480 

29 

580 

21 

500 

31 

600 

23 

520 

33 

620 

25 

540 

35 

640 

Table  3.19  Legend,  Case  (3)  and  Case  (4),  10  nm  SNR  graphs 
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S.4-S.1  Ultraviolet  comparison.  Figure  3.6  (Case  (3),  SNR  (UV), 
10  nm)  and  Figure  3.7  (Case  (4),  SNR  (UV),  10  nm)  show  the  two  graphs  compared 
in  this  section.  Table  3.20  shows  the  combined  line  orders  for  Case  (3)  and  Case  (4), 
UV  region.  The  line  spacing  comparison  is  made  using  Table  3.20,  Figure  3.6,  and 
Figure  3.7. 


ULTRAVIOLET  (UV) 

ORDER 

20 

nm 

10  nm 

5  nm 

gsma 

Top  6 

360 

360 

380 

380 

380 

5 

380 

380 

360 

360 

360 

4 

340 

400 

400 

400 

400 

3 

400 

340 

340 

340 

340 

2 

320 

320 

320 

320 

320 

320 

Bottom  1 

300 

300 

300 

300 

Missing  horizontal  lines  indicate  overlapping  wavelength  intervals 


Table  3.20  Line  Order,  Case  (3)  and  Case  (4)  Combined,  UV 


Area 

Signal  Strength 
Line  Order 
Line  Spacing 


Comparison 

Case  (3)  has  larger  SNR  values  than  Case  (4). 
Line  order  is  nearly  the  same  for  the  two  cases. 
Line  spacing  is  slightly  different. 


The  main  difference  between  Case  (3)  and  Case  (4)  in  the  Ultraviolet  region  is  in 
signal  strength.  Line  order  is  the  same  for  the  10  nm  and  5  nm  cases.  The  20  nm 
graphs  show  slight  differences  in  line  order.  These  differences  may  be  caused  by 
the  increased  bandwidth  and  increased  coarseness  of  the  spectral  reflectance.  The 
line  spacing  between  the  two  cases  is  only  slightly  different,  with  Case(4)  having  a 
number  of  lines  that  overlap.  The  Space  Object  Identification  information  for  Case 
(3)  and  Case  (4)  UV  lies  primarily  in  the  signal  strength. 


3. 4-3.2  Blue  comparison.  Figure  3.8  (Case  (3),  SNR  (B)  10  nm)  and 
Figure  3.9  (Case  (4),  SNR  (B),  10  nm)  show  the  two  graphs  compared  in  this  section. 
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SNR(1) 

SNR(3) 

SNR(5) 

^(7) 

SNR(9) 

SNR(11) 


Figure  3.6  Case  (3),  SNR  vs  Q{radm^)  (UV),  10  nm 
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Table  3.21  shows  the  combined  line  orders  for  Case  (3)  and  Case  (4)  Blue  region. 
The  line  spacing  comparison  is  made  using  Table  3.21,  Figiire  3.8,  and  Figure  3.9. 


BLUE  (B) 

ORDER 

20  nm 

10  nm 
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Top  9 
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540 
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8 
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7 
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6 
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5 
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500 
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4 
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480 

480 

3 
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460 

420 

420 

2 

420 

460 

460 

Bottom  1 

440 

440 

440 

440 

440 

Missing  horizontal  lines  indicate  overlapping  wavelength  intervals 


Table  3.21  Line  Order,  Case  (3)  and  Case  (4)  Combined,  B 


Area 

Signal  Strength 
Line  Order 
Line  Spacing 


Comparison 

Case  (3)  has  larger  SNR  values  than  Case  (4). 
Line  order  is  different. 

Line  spacing  is  different. 


Differences  exist  in  signal  strength,  line  order,  and  line  spacing  between  Case  (3) 
and  Case  (4)  in  the  B  region.  These  differences  apply  to  2dl  wavelength  bandwidths. 
Space  Object  Identification  information  may  be  stored  in  the  different  target  signal 
strengths  in  different  spectral  bands.  The  differences  between  Case  (3)  and  Case  (4) 
are  from  the  different  sizes  and  fractional  material  areas. 

S.4-S.S  Visible  comparison.  Figure  3.10  (Case  (3),  SNR  (V)  10  nm) 
and  Figure  3.11  (Case  (4),  SNR  (V),  10  nm)  show  the  two  graphs  compared  in 
this  section.  Table  3.22  shows  the  combined  line  orders  for  Case  (3)  and  Case  (4), 
Visible  region.  The  line  spacing  comparison  is  made  using  Table  3.22,  Figure  3.10, 
and  Figure  3.11. 


3-31 


Figure  3.8  Case  (3),  SNR  vs  Q{radm^)  (B),  10  nm 
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Figure  3.9  Case  (4),  SNR  vs  ft  {rad m^)  (B),  10  nm 
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VISIBLE  (V) 
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Missing  horizontal  lines  indicate  overlapping  wavelength  interv^s 


Table  3.22  Line  Order,  Case  (3)  and  Case  (4)  Combined,  V 


Area 

Signal  Strength 
Line  Order 
Line  Spacing 


Comparison 

Case  (3)  has  larger  SNR  values  than  Case  (4). 
Line  order  is  exactly  the  same  for  the  two  cases. 
Line  spacing  is  slightly  different. 


The  main  difference  between  Case  (3)  and  Caise  (4)  in  the  UV  region  is  in  signal 
strength.  The  line  spacing  between  the  two  cases  is  only  slightly  different,  with 
Case(4)  having  lines  grouped  closer  together.  The  Spaw:e  Object  Identification  infor¬ 
mation  for  Case  (3)  and  Case  (4)  Visible  region  lies  primarily  in  the  signed  strength. 

The  blue  region  offers  the  most  promise  for  discriminating  between  the  two 
targets.  The  increased  areas  of  the  solar  array,  silver  thermal  band,  and  antenna 
probably  caused  the  shift  of  the  380  nm  and  520  nm  lines.  The  380  nm  region  marks 
a  high  rate  of  change  for  the  spectral  reflectance  of  the  three  materials  (see  Figures 
F.2,  F.4,  and  F.5).  For  SOI  analysis,  looking  for  these  regions  of  dramatic  change  in 
spectral  reflectance  would  yield  the  most  hyperspectral  information. 
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SNR(19) 

SNR(21) 

SNR(23) 

SNR(25) 

SNR(27) 

SNR(29) 

^(31) 

SNR(33) 

SNR(35) 


1956.722 

1866.949 

1777.176 

1687.403 

1597.629 

1507.856 

1418.083 

1328.31 

1238.537 

1148.763 

1058.99 

969.217 

879.444 

789.67 

699.897 

610.124 

520.351 

430.578 

340.804 

251.031 

161.258 


Figure  3.10  Case  (3),  SNR  vs  U{radm^)  (V),  10  nm 


3-35 


3-36 


All  four  cases  h2Ki  enough  signal  present  for  hyperspectral  analysis  in  the  three 
wavelength  bandwidths  tested.  Differences  between  the  cases  exist  in  signal  strength, 
line  order,  and  line  spacing.  The  signal  strength,  line  order,  and  line  spacing  depend 
on  the  wavelength  and  wavelength  bandwidth  for  a  target.  Differences  in  these  areas 
indicate  Space  Object  Identification  information  may  be  stored  in  the  different  taurget 
signal  strengths  in  different  spectral  bands,  and  may  serve  as  the  foundation  for  SOI 
applications.  The  differences  existed  for  all  three  of  the  wavelength  bandwidths 
tested.  The  20  nm  bandwidth  was  consistently  different  from  the  10  nm  and  5  nm 
plots.  Even  so,  the  20  nm  bandwidth  revealed  differences  between  the  case  studies. 
For  narrow  wavelength  bandwidths,  lines  for  adjacent  regions  often  overlapped. 
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IV.  Conclusions  and  Recommendations 


4.1  Conclusions 

HyperspectraJ  analysis  is  being  proposed  as  a  means  of  enhancing  Space  Object 
Identification  capability.  This  study  developed  a  mathematical  simulation  based  on 
fundamental  radiometry  principles  to  evaluate  the  feasibility  of  hyperspectral  anal¬ 
ysis  from  a  Signal  to  Noise  perspective.  For  the  SNR  evaluation,  the  target  signal, 
background  and  detector  noise,  and  integration  time  had  to  be  characterized.  The 
target  signal  was  characterized  by  the  photon  flux,  spectral  reflectivity,  and  a  hy- 
pernormalized  parameter  that  grouped  all  remaining  variables  into  a  single  variable. 
The  background  noise  was  assumed  small  in  comparison  to  the  target  signal,  and 
detector  noise  was  modeled  from  a  CCD  detector.  The  integration  time  was  bounded 
using  a  best  and  worst  case.  The  best  case  was  applied  to  the  radiometry  model. 
The  worst  case  was  determined  by  selecting  a  minimum  SNR  threshold  and  solving 
for  the  integration  time.  Using  these  parameters,  the  single  pixel  SNR  equation 
for  hyperspectrometry  and  SNRp^j.  equation  for  hyperspectr2d  imaging  were  devel¬ 
oped.  The  count  rate  from  the  detector,  the  single  pixel  SNR,  and  SNRp^j.  are 
outputs  of  the  radiometry  model.  The  radiometry  model  was  used  for  four  case 
studies.  Two  cases  involved  hyperspectred  imaging  and  two  involved  hyperspectrom¬ 
etry.  The  procedure  for  running  a  case  study  involved  specifying  a  configuration, 
interpolating  the  spectral  reflectances,  and  running  the  radiometry  model.  For  the 
hyperspectral  imaging  Ceise  studies,  the  modeled  target  was  a  rocket  body,  made  of 
aluminum  in  Case  (1),  and  painted  white  in  Case  (2).  For  the  hyperspectrometry 
case  studies  Case  (3)  and  Case  (4),  the  targets  were  two  communication  satellites, 
having  the  same  shape  but  different  sizes  and  fractional  material  areas.  The  case 
studies  revealed  the  following  main  conclusions: 

1.  Enough  signal  is  present  for  hyperspectral  imaging  and  hyperspectrometry  for 

wavelength  bandwidths  of  A  A=20  nm,10  nm,  and  5  nm. 
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2.  The  results  for  Case  (1)  and  Case  (2),  and  for  Case  (3)  and  Case  (4),  show 
differences  in  sign2d  level  as  a  function  of  wavelength  and  bandwidth  which 
may  serve  as  a  basis  for  Space  Object  Identification.  For  Case  (1)  and  Case 
(2),  differences  can  be  attributed  to  the  different  spectral  reflectances  between 
an  aluminum  and  white  object.  For  Case  (3)  and  Case  (4),  the  difference  can 
be  attributed  to  the  different  sizes  zmd  fractional  material  areas. 

These  conclusions  prove  the  feasibility  of  hyperspectral  measurements  of  space  ob¬ 
jects.  A  number  of  corollary  issues  regarding  hyperspectral  analysis  were  observed 
in  the  course  of  this  analysis: 

1.  Two  or  more  lines  often  overlapped  in  the  count  rate  and  SNR  graphs.  For  the 
visible  region  in  Case  (3)  and  Case  (4),  three  wavelength  intervals  (600  nm, 
620  nm,  and  640  nm)  always  overlapped.  This  problem  is  compounded  when 
two  regions  are  adjacent,  especially  at  AA  =  5  nm. 

2.  All  the  wavelength  bandwidths  showed  differences  between  the  targets.  If 
a  bandwidth  of  AA=:  20  nm  reveals  significant  SOI  characteristics,  making 
the  wavelength  bandwidth  more  narrow  is  not  necessary.  For  data  loading 
considerations,  a  AA=5  nm  requires  at  least  four  times  the  data  processing 
capability  as  a  AA=20  nm. 

J^.2  Recommendations  for  Further  Research 

This  study  demonstrated  enough  target  signal  is  present  to  perform  hyperspec¬ 
tral  analysis.  The  following  areas  of  research  are  offered  ais  possibilities  for  further 
investigation: 

1.  Given  a  count  rate  from  a  detector,  can  the  target’s  material  composition  be 
determined? 

2.  Develop  a  hyperspectral  imaging  simulation  that  links  a  target’s  spectral  re¬ 
flectance  to  spatial  coordinates. 
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3.  For  a  ground  based  sensor,  investigate  atmospheric  effects  on  hyperspectral 
analysis. 

4.  What  is  the  maximum  wavelength  bandwidth  that  would  provide  SOI  infor¬ 
mation,  and  keep  the  data  load  at  a  minimum?  Are  there  regions  of  the 
electromagnetic  spectrum  that  reveal  more  SOI  information  than  others?  For 
exjunple,  in  Case  (3)  and  Case  (4),  the  blue  region  revealed  the  greatest  differ¬ 
ences  between  the  two  targets. 

In  summary,  enough  signal  is  present  for  hyperspectral  analysis  of  space  objects.  The 
case  studies  in  this  analysis  show  differences  exist  between  targets  which  might  be 
exploited  for  SOI  purposes.  Hyperspectral  analysis  warrants  further  investigation. 
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Appendix  A.  Mathcad  Models 


A.l  Model  Configuration 

Before  running  a  case  study,  two  things  must  be  determined;  the  ability  of 
the  sensor  to  image  the  target  and  the  reference  hypemormalized  parameter.  The 
Mathcad  configmation  model  (page  A-3)  calculates  the  target  angular  subtense  and 
angular  subtense  of  the  diffraction  limited  spot.  When  <  ^4>sah  sensor 

can  image  the  target.  Using  the  target  and  sensor  parameters,  the  hypemormalized 
parameter  fl  is  calculated  for  normal  incidence.  This  value  becomes  the  hypernor- 
malized  reference  parameter,  fire/- 

A. 2  Spectral  Reflectivity  Interpolation 

The  spectral  reflectances  for  a  material  are  compiled  in  the  Matter.dat  Surface 
Properties  Data  Compilation  (14).  This  compilation  assigns  a  Surface  Material  Code 
(SMC)  to  a  material  and  gives  the  spectral  reflectances  versus  wavelength  in  incre¬ 
ments  of  10  nm.  The  Mathcad  interpolation  code  (pages  A-4  and  A-6)  generates 
spectral  reflectances  at  20  nm  and  5  nm  bandwidths.  The  code  also  produces  plots 
of  the  spectral  reflectances  for  each  bandwidth  to  verify  the  interpolated  values  give 
a  spectral  reflectance  curve  close  to  the  original  spectral  reflectance  curve  in  the 
Matter.dat  Surface  Properties  Data  Compilation, 

A.S  Radiometry  Model 

The  Mathcad  radiometry  code  (page  A-7)  generates  the  count  rate,  SNR,  and 
SNRp^j.  plots.  The  model  requires  the  maximum  and  minimum  wavelengths  and 
wavelength  bandwidth  (A  A)  to  calculate  the  photon  flux  P  using  Planck’s  formula 
and  the  inverse  square  law.  A  counter  (i)  is  used  as  a  reference  pointer  throughout 
the  radiometry  model,  and  corresponds  to  a  wavelength  interval.  The  last  page  of  this 
appendix  (page  A-8)  contains  a  counter  to  wavelength  conversion  for  20  nm,  10  nm 
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and  5  nm  wavelength  bandwidths.  For  example,  if  A  A  =  10  nm,  an  t  =  23  is  a  10  nm 
bandwidth  centered  on  520  nm.  Spectral  relectance  data  files  are  specific  to  material 
and  wavelength.  The  data  file  associated  with  po2  for  a  20  nm  A  A  is  smc02xx.dat, 
XX  for  20  nm  (x  for  10  nm,  and  v  for  5  nm).  The  hypernormalized  parameter  is 
varied  by  an  order  of  magnitude  above  and  below  the  reference  hypernormalized 
parameter’s  exponent.  The  remaining  input  parameters  are  the  integration  time  td, 
and  the  detector  read  noise  arn-  Mathcad  calculates  the  count  rate,  the  overall  SNR 
and  the  SNRp^^  (if  applicable),  and  outputs  the  results  to  graphs. 
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Configuntion 

Target’  Rocketbody,6mkng,2mmdUmeier.  RangeSOOkm. 
Case  1 A1  and  case  2  !)^te  paint 
Sensor  Diain^erprinuuy  mirror  «  1.6  m. 


Target  parameters: 

Sensor  Parameters 

L  j  6*m  Length 

X  :=  900‘iim  Max  wavdengdi 

r  :=  800>km  Range 

D:=liS‘m  Sensor  optic  diameter 

L 

^max 

d4^„:=l.22. 

^mm  p 

^min^^sat 

Can  image  target 

C(xiq)ute  hypemonnalized  parameter  diat  conesponds  to  above  object’ 
Hypemonnalized  parameter 

c  I^V 

5  R  a*l  yl  Area  of  sensor  opbc 

rs800*km  Range  to  target 

6  j  :=  0  Bound  hypemonnalized  parameter  at  normal  incideoce 

62-0 

X  o  -  -98  Optical  transmissitm 

S  j  :=  12-mP  Area  of  target 

Ti  :=  .6  Average  quantum  efficiency  detector 
*R 

“  ref  -  — «»(0 1)  •«»(e  2)  ’X  ^  S 


Q 


ref 


!  2.217*  10  *m^  Reference  hypemonnalized  parameter 


-3 


F 


Liteipolitioii 

This  file  tikes  data  from  a  .dat  fik  built  from  the  Matter.dat  Swfrice  Properties  Data 
Coaqtilatkm,  and  mterpoiates  to  build  q)ectral  reflectance  vectors  of  20  nm 
and  S  nm  wavelength  bandwidths. 

i:=  1..69 

vxj  :=  i  i  vx.  vy.  x  G)unters  and  variables  used  to  do  the  inierpdation 


vy.  R£AD(smcl8x) 
X  :=  1 ..  69 


Reads  in  data  from  the  riat  file.  To  link  variaUe  name 
to  filename,  use  "Associate  Filename..."  under  File. 


69 

10  nm  (Original  wavelength  bandwidth) 


Verify  spectral  reflectance  grqrh  is  q)pioximately  equivalent  to  original 
the  gnq)h  in  the  Matter  jdat  Surface  Properties  Data  Compilation. 
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20  nm  interpoUtioiL  The  x-  is  multiplied  by  the  increment  of  2,  — 
where  20  is  the  desired  dX  and  10  is  the  dX  read  in  above, 
i  :=  I ..  35 


WRITE(smcl8xx)  :=  f 

Writes  to  INTRP^smc??xx.dat 
Change  name  to  match  SMC#. 


f  :s  linteip(vx,vy,x) 
length(f)  s36 


20  nm  Interpolation 

Verify  spectral  reflectance  graph  is  apimndmately  equivalent  to  the 
(xiginal  gnq)h  in  the  Matter  Surface  Properties  Data  Compilation. 
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is  mult^lied  by  the  increment  of 
and  10  is  the  dX  read  in  abow. 


S  nm  intopolation.  The 
where  S  is  tlw  desired  dX 

i  :=  1 ..  141 

.  5 
X.  :=  i~-“ 

'  10 

- » 

f  :=  linteip(vx,vy.x) 

length(f)  s  142 


WRlTE(smcl8v)  :=  f 

Writes  to  INTRP^smc??v.dat 
Change  namg  to  match  SMC#. 


*.,vx.  70.5 

5  nm  Interpolation 


Verify  spectral  reflectance  graph  is  iq;>proximately  equivalent  to  the 
original  gr^h  in  the  Matter  jia  Surface  Properties  Data  Compilation. 
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Rndtoroetry 

X  jb  •*  290'mn  X  jjb  •=  910-nm  Wavelength  lower  and  ui^ier  bound 

dX  20-nm  Wavelength  bandwidth 


„  ^ub'^lb 

N  :* - 

AX 

i:=l..N 
b.  :=i  AX  +  XiK 


a.  :*  bj  -  AX 


N  «  31  Number  oi  wavelength  bands 
Counts 

Upper  limit  inciementer  for  integral 
Lower  limit  incrementer  for  integral 


2-iC’h‘C* 


X’.  . 


Planck's  formula 


r*>i 

My,(X) 

- XdX 

p  /”^sunj 

i  •'  1 

h-c-jc 

\  sun/ 

«i 

Photon  flux 


P  02.  READ(stnc02xx)  Spectral  reflectivity 
Q  :=  l(r«.m’.i(r“-»-m*..  l(r‘®mi*  Rangeof  Q 


ix(i)  q  Q-Pi  Po2. 

t^  8-sec  bitegration  time 

SNR(i)  :=  --  ?- 

iT(i)  io 

- -t-  — 

N  q  q 


COUNT  RATE  CALCULATION 
Deiector 

<5  _  :=  12  “«se 

ni 

SNR  CALCULATION 


a.i-b. 


SNRpix(i):= 


Mean  wavelength  z 0,.0l ..  li)  Angles  parameter 
*^i‘P02.’*d 

=  :  SNR  ^  CALCULATION 
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20  an  biBdwkUl 
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11 
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980 
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440 
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28 
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9 
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19 

660 

29 
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10 

480 

20 

680 

30 

880 

10  am  bttidwidth 


1 

300 

11 

400 

21 

500 

31 

600 

41 

700 

51 

800 

61 

900 

2 

310 

12 

410 

22 

510 

32 

610 

42 

710 

52 

810 

62 

910 

3 

320 

13 

420 

23 

520 

33 

620 

43 

720 

53 

820 

63 

920 

4 

330 

14 

430 

24 

530 

34 

630 

44 

730 

54 

830 

64 

930 

5 

340 

15 

440 

25 

540 

35 

640 

45 

740 

55 

840 

65 

940 

6 

350 

16 

450 

26 

550 

36 

650 

46 

750 

56 

850 

66 

950 

7 
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17 

460 

27 

560 

37 

660 

47 

760 

57 

860 

67 

960 

8 

370 

18 

470 

28 

570 

38 

670 

48 

770 

58 

870 

68 

970 

9 

380 

19 

480 

29 

580 

39 

680 

49 

780 

59 

880 

69 

980 

10 

390 

20 

490 

30 

590 

40 

690 

50 

790 

60 
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70 

990 

5  nm  bandwidth 


1  300 

11  350 

21  400 

31  450 

41  500 

51  550 

61  600 

2  305 

12  355 

22  405 

32  455 

42  505 

52  555 

62  60S 

3  310 

13  360 

23  410 

33  460 

43  510 

53  560 

63  610 

4  315 

14  365 

24  415 

34  465 

44  515 

54  565 

64  615 

5  320 

15  370 

25  420 

35  470 

45  520 

55  570 

65  620 

6  325 

16  375 

26  425 

36  475 

46  525 

56  575 

66  625 

7  330 

17  380 

27  430 

37  480 

47  530 

57  580 

67  630 

8  335 

18  385 

28  435 

38  485 

48  535 

58  585 

68  635 

9  340 

19  390 

29  440 

39  490 

49  540 

59  590 

69  640 

10  345 

20  395 

30  445 

40  495 

50  545 

60  595 

70  645 

71  650 

81  700 

91  750 

101  800 

111  850 

121  900 

131  950 

72  655 

82  70S 

92  755 

102  80S 

112  855 

122  905 

132  955 

73  660 

83  710 

93  760 

103  810 

113  860 

123  910 

133  960 

74  665 

84  715 

94  765 

104  815 

114  865 

124  915 

134  965 

75  670 

85  720 

95  770 

105  820 

115  870 

125  920 

135  970 

76  675 

86  725 

96  775 

106  825 

116  875 

126  925 

136  975 

77  680 

87  730 

97  780 

107  830 

117  880 

127  930 

137  980 

78  685 

88  735 

98  785 

108  835 

118  885 

128  935 

138  985 

79  690 

89  740 

99  790 

109  840 

119  890 

129  940 

139  990 

80  695 

90  745 

100  795 

110  845 

120  895 

130  945 

140  995 

A-8 


Appendix  B.  Case  (1):  20  nm,  JO  nm  and  5  nm  Graphs 

B.l  Legend 

The  legend  for  Case  (1)  graphs  is  given  in  Table  B.l. 
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Table  B.l  Legend,  Case  (1) 


B.2  20  nm  graphs 

Figure  B.l  shows  the  target  signal  count  rate  ij/q  vs  Q  for  Case  (1),  AA  = 
20 nm.  The  vertical  dashed  line  is  the  reference  hypernormalized  parameter,  fire/, 
computed  for  the  target. 

Figure  B.2  shows  the  SNR  vs  Q  for  Case  (1),  A  A  =  20  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  fire/)  computed  for  the  target. 

Figure  B.3  shows  the  SNRp^j.  vs  z  for  Case  (1),  AA  =  20 nm. 
b.S  10  nm  graphs 

Figure  B.4  shows  the  target  signal  count  rate  ix/q  vs  fi  for  Case  (1),  A  A  = 
10  nm.  The  vertical  dashed  line  is  the  reference  hypernormalized  parameter,  flrefi 
computed  for  the  target. 

Figure  B.5  shows  the  SNR  vs  ft  for  Case  (1),  A  A  =  10  nm.  The  verticad  dashed 
line  is  the  reference  hypernormalized  parameter,  ftre/,  computed  for  the  target. 

Figure  B.6  shows  the  SNR—j.  vs  z  for  Case  (1),  AA  =  10 nm. 
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B.4  5  nm  graphs 


Figure  B.7  shows  the  target  signal  count  rate  ij/q  vs  Cl  for  Case  (1),  AA  = 
5nm.  The  vertical  dashed  line  is  the  reference  hypernormalized  p2irameter,  fire/, 
computed  for  the  target. 

Figure  B.8  shows  the  SNR  vs  Cl  for  Case  (1),  A  A  =  5nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

Figure  B.9  shows  the  SNR—^.  vs  z  for  Case  (1),  A  A  =  5nm. 
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Case  (1),  SNR  vs  S}(radm^),  20  nm 
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Figure  B.3  Case  (1),  SNRjj^j.  vs  z,  20  nm 
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Figure  B.5  Case  (1),  SNR  vs  fi  {radm'^),  10  nm 
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Figure  B.6  Case  (1),  SNRj^j.  vs  z,  10  nm 
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Figure  B.8  Case  (1),  SNR  vs  (rodm*),  5  nm 
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Appendix  C.  Case  (2):  20  nm,  10  nm  and  5  nm  Graphs 

C.l  Legend 

The  legend  for  Case  (2)  graphs  is  given  in  Table  C.l. 
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Table  C.l  Legend,  Case  (2) 


C.2  20  nm  graphs 

Figure  C.l  shows  the  target  signal  count  rate  ij/q  vs  0,  for  Case  (2),  AA  = 
20 nm.  The  vertical  dashed  line  is  the  reference  hypernormalized  parameter,  fire/, 
computed  for  the  target. 

Figure  C.2  shows  the  SNR  vs  fl  for  Case  (2),  A  A  =  20  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

Figure  C.3  shows  the  SNRp^j.  vs  z  for  Case  (2),  A  A  =  20  nm. 

C.3  10  nm  graphs 

Figure  C.4  shows  the  target  signal  count  rate  ir/q  vs  fl  for  Case  (2),  AA  = 
10 nm.  The  vertical  dashed  line  is  the  reference  hypernormalized  parameter,  fire/, 
computed  for  the  target. 

Figure  C.5  shows  tne  SNR  vs  fl  for  Case  (2),  A  A  =  10  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  p£irameter,  fire/,  computed  for  the  target. 

Figure  C.6  shows  the  SNR~^  vs  z  for  Case  (2),  AA  =  10 nm. 
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C.4  5  nm  graphs 


Figure  C.7  shows  the  target  signal  count  rate  ij/q  vs  Q  for  Case  (2),  AA  = 
5nm.  The  vertical  dashed  line  is  the  reference  hypemormalized  parameter,  ilref, 
computed  for  the  target. 

Figure  C.8  shows  the  SNR  vs  fl  for  Case  (2),  A  A  =  5uxn.  The  vertical  dashed 
line  is  the  reference  hypemormalized  parameter,  fire/i  computed  for  the  target. 

Figure  C.9  shows  the  SNRj^^  vs  z  for  Case  (2),  AA  =  5nm. 
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Figure  C.3  Case  (2),  SNR-jj.  vs  z,  20  nm 
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Figure  C.6  Case  (2),  SNE^x  vs  z,  10  nm 
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Figure  C.7  Case  (2).  Count  Rate  ^  (1/s)  vs  U{radm^),  5  nm 
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Appendix  D.  Case  (3):  20  nm,  10  nm  and  5  nm  Graphs 
D.  1  Legend 

Table  D.l  shows  the  legend  for  Case  (3)  graphs. 
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Table  D.l  Legend,  Case  (3) 
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D.2  20  nm  graphs 

D.2.1  Ultraviolet  (UV),  300  nm  -  4OO  nm.  Figure  D.l  shows  the  target 
signal  count  rate  ij/q  vs  Q  for  Case  (3),  AA  =  20 nm.  The  vertical  dashed  line  is 
the  reference  hypernormalized  parameter,  Urefi  computed  for  the  target. 

Figure  D.2  shows  the  SNR  vs  fi  for  Case  (3),  A  A  =  20  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  Dre/>  computed  for  the  target. 

D.2. 2  Blue  (B),  380  nm  -  54O  nm.  Figure  D.3  shows  the  target  signal 
count  rate  ij/q  vs  fl  for  Case  (3),  A  A  =  20  nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  Dr*/,  computed  for  the  target. 

Figure  D.4  shows  the  SNR  vs  Q  for  Case  (3),  A  A  =  20  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

D.2. 3  Visible  (V),  4^0  nm  -  650  nm.  Figure  D.5  shows  the  target  signal 
count  rate  ij/q  vs  fl  for  Case  (3),  AA  =  20 nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  flre/>  computed  for  the  target. 

Figure  D.6  shows  the  SNR  vs  fl  for  Cawe  (3),  A  A  =  20  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

D.S  10  nm  graphs 

D.3.1  Ultraviolet  (UV),  300  nm  -  4OO  nm.  Figure  D.7  shows  the  target 
signal  count  rate  ij/q  vs  fl  for  Case  (3),  AA  =  10 nm.  The  vertical  dashed  line  is 
the  reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

Figure  D.S  shows  the  SNR  vs  fl  for  Caise  (3),  A  A  =  10  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 
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D.S.2  Blue  (B),  380  nm  -  54O  nm.  Figure  D.9  shows  the  taurget  signed 
count  rate  ij/q  vs  for  Case  (3),  AA  =  10 nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

Figure  D.IO  shows  the  SNR  vs  f)  for  Case  (3),  AA  =  10 nm.  The  vertical 
daished  line  is  the  reference  hypernormalized  parameter,  fire/,  computed  for  the 
target. 

D.S.S  Visible  (V),  480  nm  -  650  nm.  Figure  D.ll  shows  the  target  signal 
count  rate  ij/q  vs  fi  for  Case  (3),  AA  =  10 nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

Figure  D.12  shows  the  SNR  vs  fi  for  Case  (3),  AA  =  10 nm.  The  vertical 
dashed  line  is  the  reference  hypernormalized  parameter,  fire/,  computed  for  the 
target. 

D.4  5  nm  graphs 

D.4‘1  Ultraviolet  (UV),  300  nm  -  4OO  nm.  Figure  D,13  shows  the  target 
signal  count  rate  ij/q  vs  fi  for  Case  (3),  AA  =  5nm.  The  vertical  dashed  line  is 
the  reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

Figure  D.14  shows  the  SNR  vs  fi  for  Case  (3),  A  A  =  5  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

D.4-2  Blue  (B),  380  nm  -  54O  nm.  Figure  D.15  shows  the  target  signal 
count  rate  ij/q  vs  fi  for  Case  (3),  AA  =  5nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

Figure  D. 16  shows  the  SNR  vs  fi  for  Case  (3),  A  A  =  5nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 
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D.4-S  Visible  (V),  48O  nm  -  650  nm.  Figure  D.17  shows  the  target  signal 
count  rate  ij/q  vs  fi  for  Case  (3),  AA  =  5nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  flre/i  computed  for  the  target. 

FigureD.lSshowstheSNRvsnfor  Case(3),  AA  =  5nm.  The  verticad  dashed 
line  is  the  reference  hypernormalized  parameter,  fir*/,  computed  for  the  target. 
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Figure  D.2  Case  (3),  SNR  vs  U{radm'^)  (UV),  20  nm 
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Figure  D.3  Case  (3),  Count  Rate  ^  (1/s)  vs  ^(radm^)  (B),  20  nm 
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Figure  D.5  Case  (3),  Count  Rate  ^  (l/s)  vs  n{radm^)  (V),  20  nm 
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Figure  D.6  Case  (3),  SNR  vs  (radm^)  (V),  20  nm 
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Appendix  E.  Case  (4)'-  20  nm,  10  nm  and  5  nm  Graphs 

E.l  Legend 

Table  E.l  shows  the  legend  for  Case  (4). 


Table  E.l  Legend,  Case  (4) 
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E.2  20  nm  graphs 

E.2.1  Ultraviolet  (UV),  $00  nm  -  4OO  nm.  Figure  E.l  shows  the  target 
signal  count  rate  ij/q  vs  fi  for  Case  (4),  AA  =  20 nm.  The  vertical  dashed  line  is 
the  reference  hypernormalized  parameter,  computed  for  the  target. 

Figure  E.2  shows  the  SNR  vs  n  for  Case  (4),  A  A  =  20  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  fire/)  computed  for  the  target. 

E.2. 2  Blue  (B),  380  nm  -  54O  nm.  Figure  E.3  shows  the  target  signal 
count  rate  ijlq  vs  for  Case  (4),  AA  =  20 nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  n,e/>  computed  for  the  target. 

Figure  E.4  shows  the  SNR  vs  fl  for  Case  (4),  A  A  =  20  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

E.2. 3  Visible  (V),  48O  nm  -  650  nm.  Figure  E.5  shows  the  target  signal 
count  rate  ijlq  vs  fl  for  Case  (4),  AA  =  20 nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  flre/j  computed  for  the  target. 

Figure  E.6  shows  the  SNR  vs  H  for  Case  (4),  A  A  =  20  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  pcirameter,  fire/?  computed  for  the  target. 

E.3  10  nm  graphs 

E.3.1  Ultraviolet  (UV),  300  nm  -  4OO  nm.  Figure  E.7  shows  the  target 
signal  count  rate  ijfq  vs  U  for  Case  (4),  AA  =  10 nm.  The  vertical  dashed  line  is 
the  reference  hypernormalized  parameter,  tlref,  computed  for  the  target. 

Figure  E.8  shows  the  SNR  vs  17  for  Case  (4),  A  A  =  10  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  ilre/f  computed  for  the  target. 
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E.S.2  Blue  (B),  380  nm  -  54O  nm.  Figure  E.9  shows  the  target  signal 
count  rate  ij/q  vs  ft  for  Case  (4),  AA  =  10 nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  ft,*/,  computed  for  the  target. 

Figure  E.IO  shows  the  SNR  vs  ft  for  Case  (4),  AA  =  10 nm.  The  vertical 
dashed  line  is  the  reference  hypernormalized  parameter,  ftre/,  computed  for  the 
target. 

E.S.S  Visible  fVJ,  4^0  nm  -  650  nm.  Figure  E.ll  shows  the  target  signal 
count  rate  ir/q  vs  ft  for  Case  (4),  AA  =  10 nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  ftre/,  computed  for  the  target. 

Figure  E.12  shows  the  SNR  vs  ft  for  Case  (4),  AA  =  10 nm.  The  vertical 
dashed  line  is  the  reference  hypernormalized  parameter,  ftre/,  computed  for  the 
target. 

E.4  5  nm  graphs 

E.4'1  Ultraviolet  (UV),  300  nm  -  4OO  nm.  Figure  E.13  shows  the  target 
signal  count  rate  ij/q  vs  ft  for  Case  (4),  AA  =  5nm.  The  vertical  dashed  line  is 
the  reference  hypernormalized  parameter,  ftre/,  computed  for  the  target. 

Figure  E.14  shows  the  SNR  vs  ft  for  Case  (4),  A  A  =  5  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  ftre/,  computed  for  the  target. 

E.4-2  Blue  (B),  380  nm  ~  54O  nm.  Figure  E.15  shows  the  target  signal 
count  rate  ix/q  vs  ft  for  Case  (4),  AA  =  5nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  ftre/,  computed  for  the  target. 

Figure  E.16  shows  the  SNR  vs  ft  for  Case  (4),  A  A  =  5  nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  ftre/,  computed  for  the  target. 
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E.4‘S  Visible  (V),  480  nm  -  650  nm.  Figure  E.17  shows  the  taxget  signal 
count  rate  ij/q  vs  fl  for  Case  (4),  AA  =  5nm.  The  vertical  dashed  line  is  the 
reference  hypernormalized  parameter,  fire/,  computed  for  the  target. 

Figure  E. 18  shows  the  SNR  vs  f)  for  Case  (4),  A  A  =:  5nm.  The  vertical  dashed 
line  is  the  reference  hypernormalized  parameter,  computed  for  the  target. 
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Figure  E.8  Case  (4),  SNR  vs  O(rarfm^)  (UV),  10  nm 
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Figure  E.ll  Case  (4),  Count  Rate  ^  (1/s)  vs  (V),  10  nm 
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Figure  E.15  Case  (4),  Count  Rate  ^  {l/s)  vs  U{radm^)  (B),  5  nm 
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Figure  E.16  Case  (4),  SNR  vs  Q,{radm})  (B),  5  nm 
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Figure  E,18  Case  (4),  SNR  vs  Q{radm^)  (V),  5  nm 
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Appendix  F.  Spectral  Reflectance  of  Various  Materials 

Figure  F.l  shows  the  spectraJ  reflectance  versus  wavelength  for  SMC-02,  Alu¬ 
minum  Alloy  Type  2024-T3  (Mill  Finish). 

Figure  F.2  shows  the  spectral  reflectance  versus  wavelength  for  SMC-18,  Silver. 

Figure  F.3  shows  the  spectral  reflectance  versus  wavelength  for  SMC-20,  Paint, 
Chemglaze  A276  (White). 

Figure  F.4  shows  the  spectral  reflectance  versus  wavelength  for  SMC-23,  Solar 
Cell,  Silicon. 

Figure  F.5  shows  the  spectral  reflectance  versus  wavelength  for  SMC-30,  Mylar, 
Aluminized  (Aluminum  Side  Out), 
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Figure  F.2  Spectral  Reflectance  (%),  SMC-18:  Silver 
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Figure  F.5  Spectral  Reflectance  (%),  SMC*30:  Mylar,  Aluminized  (Aluminum  Side 
Out) 
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